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FOREWORD 


This  report,  and  the  associated  Secret  Supplement, 
documents  the  third  phase  of  the  effort  performed  for  the  Naval 
Sea  Systems  Command  Code  63R-11  under  Contract  N00024-79-C- 
6405,  and  covers  the  period  from  October  1979  to  October  1980. 
This  study  has  as  its  goal  the  examination  of  adaptive  filters  to 
passive  sonar  bearing  tracking.  This  has  been  accomplished 
via  statistical  modeling  and  analysis,  with  verification  by  simu¬ 
lation  on  computer  generated  data  and  data  recorded  at  sea. 
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1.0  INTRODUCTION 

^This  is  the  Final  report  on  phase  three  of  the  Adaptive  Tracking  System  Study, 
which  considers  the  application  of  adaptive  filtering  to  passive  sonar  bearing  esti¬ 
mation.  The  study  was  initiated  because  the  adaptive  filter  structure  seemed  well 
suited  to  bearing  estimation  problems  involving  unknown  Input  statistics  and  dy¬ 
namic  targets.  In  particular,  (a)  an  adaptive  filter  does  not  require  a  priori  power 
spectral  information  on  signal  and  noise  fields,  (b)  because  it  is  iterative,  it  can 
track  non-stationary  Inputs  while  preserving  its  underlying  minimum  mean  square 
optimality  criterion,  and  (c)  since  all  the  correlation  information  between  adaptive 
filter  inputs  is  contained  in  the  filter  weights,  the  potential  exists  to  perform  both 
broadband  and  narrowband  tracking  simultaneously. 

The  first  phase  of  this  study  ^  addressed  the  feasibility  of  such  a  tracker 
using  the  Least  Mean  Square  (LMS)  adaptive  filter  algorithm,  and  developed  the 
tracker  structure  of  Figure  1.  to  this  structure,  the  two  beamfoimed  split  array 
outputs  are  provided  as  the  two  inputs  to  an  LMS  adaptive  filter  configured  as  a 
canceller.  The  output  of  one  half  array,  processed  by  a  non-recursive  adaptive 
filter,  is  subtracted  from  the  other  half  array  to  yield  an  error  signal.  The  error 
signal  is  then  used  to  recursively  update  the  adaptive  filter  weights,  or  impulse 
response,  such  that  the  mean  square  error  is  minimized.  If  x(nTs)  and  d(nTs)  are 
the  two  inputs,  with  T  the  sample  rate  and  n  the  time  index,  the  adaptive  filter 
stores  the  data  vector 

X(n)  -  [x [nTs],  x  [(n-l)Ts],.  ..x  [(n-M+l)Ts]]T  (1) 

T 

where  X  denotes  X  transpose.  The  output  of  the  adaptive  filter  is 
y(n)  *  WT(n)X(n) 

where  W(n)  is  a  vector  of  M  weights  at  time  nT_.  The  error  output  is  therefore 

mam  9 


(2) 


Figure  1.  Adaptive  Tracker  Structure 


Figure  1.  Adaptive  Tracker  Structure 
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€(n)  *  d(n)  -  WT(n)X(n). 

The  LMS  algorithm  updates  the  weight  vector  on  each  iteration  as 

W(n+1)  =  W(n)  +  u  €(n)X(n)  (4) 

where  n  is  a  weight  update  coefficient.  This  parameter  controls  the  rate  of  con¬ 
vergence,  algorithm  noise,  and,  ultimately,  the  stability  of  the  algorithm. 

In  converging  to  minimize  the  mean  square  error  between  half  array  outputs, 
the  adaptive  filter  must  incorporate  any  time  delay  (or  for  narrowband  sources, 
phase  shift)  between  its  inputs  in  the  weight  vector.  For  broadband  inputs,  the 
weights  have  the  same  shape  as  the  signal  autocorrelation  function  with  the  peak 
located  at  the  delay  between  the  phase  centers.  Hence,  the  tracker  must  determine 
the  location  of  the  peak  o f  the  weights,  using  Interpolation  between  discrete  taps  if 
necessary.  For  sinusoidal  inputs,  the  weights  converge  to  a  sinusoid  with  phase 
equal  to  the  phase  shift  between  array  halves.  The  tracker  then  determines  the 
phase  shift  by  means  of  spectrum  analysis  of  the  weight  vector.  Either  the  time  de¬ 
lay  or  phase  shift  is  easily  converted  to  a  bearing  estimate. 

The  first  phase  of  the  study  also  provided  many  of  the  analytical  techniques 
required  for  analysis  of  the  adaptive  tracker.  In  particular,  it  provided  the  fre¬ 
quency  domain  model  for  the  adaptive  tracker  shown  in  Figure  2.  In  the  frequency 
domain  version,  the  two  inputs  are  Fast  Fourier  Transformed  (FFT'd)  and  a  single 
tap,  complex  LMS  adaptive  filter  performed  Independently  to  the  time  domain  filter, 
both  in  transient  and  steady  state  operation,  under  certain  conditions  on  the  filter 
design  parameters.  The  primary  advantage  of  the  frequency  domain  model  is  that 
both  the  mean  and  variance  of  the  weights  can  be  predicted.  By  careful  selection  of 
filter  parameters,  this  model  can  then  be  used  to  predict  the  variance  of  the  time 
domain  weights,  which  cannot  be  obtained  by  analysis  in  the  time  domain. 
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Given  the  statistics  of  the  time  and  frequency  domain  weights,  it  was  possible 
to  obtain  predictions  of  the  bearing  estimate  statistics  for  broadband  inputs.  A  con¬ 


tinuous  broadband  adaptive  tracker  was  shown  to  perform  within  0.5  dB  of  the 
Cramer-Rao  Lower  Bound  on  the  variance  of  all  unbiased  estimators.  Performance 
results  were  also  derived  for  a  practical  realization  of  the  tracker,  operating  at  a 
reasonable  sampling  rate  with  interpolation  between  samples  providing  fine  delay 
(hence  bearing)  resolution.  These  results  were  shown  to  be  comparable  to  a  con¬ 
ventional  tracker.  The  first  phase  of  the  study  also  developed  the  mean  weight  of 
the  adaptive  filter  with  dynamic  broadband  and  narrowband  inputs.  The  results  were 
applied  in  Phase  Two  in  the  analysis  of  the  tracker  performance  for  dynamic  targets. 
Finally,  the  first  phase  developed  extensive  simulation  programs  for  the  tracker 
structure,  which  were  used  to  validate  the  analyses,  and  provide  an  initial  look  at 
dynamic  tracking  performance. 

The  second  phase  of  the  study  considered  in  detail  the  potential  advantages 
of  the  adaptive  tracker  with  combined  broadband  and  narrowband  inputs,  unknown 
a  priori  input  spectra,  and  dynamic  targets.  The  bearing  estimate  variance  was 
determined  analytically  for  narrowband  targets,  combined  narrowband  and  broad¬ 
band  targets,  and  narrowband  targets  in  the  presense  of  a  broadband  Interference. 
For  a  hi signal  to  noise  ratio  narrowband  target,  the  adaptive  tracker  was  shown 
to  perform  within  0. 5  dB  of  the  Cramer-Rao  lower  bound  on  the  variance  of  any 
unbiased  estimator  operating  on  the  same  observations.  It  was  also  shown  that  in 
the  presense  of  interference,  the  bias  of  the  estimate  dominates  the  bearing  error 
for  low  SNR  targets. 

The  adaptive  tracker  was  compared  to  a  conventional  split  bearing  tracker  when 
the  a  priori  assumptions  on  the  input  spectra  were  incorrect.  The  adaptive  tracker 
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was  shown  to  be  significantly  more  tolerant  of  changes  in  the  input  spectra  than 
conventional  trackers.  Whether  this  advantage  can  be  realized  in  actual  practice 
will  depend  upon  the  degree  to  which  input  spectra  vary  in  actual  sonar 
environments. 

An  analysis  comparing  the  adaptive  tracker  and  a  conventional  split  bearing 
tracker  with  moving  targets  was  also  performed.  For  a  specific  input  correlation 
function,  the  adaptive  tracker  was  shown  to  be  about  1. 6  dB  bettern  than  the  con¬ 
ventional  tracker. 

The  third  phase  of  the  study  concentrated  in  three  areas.  It  continued  the 
analytical  evaluation  of  the  tracker  to  include  such  effects  as  multipath,  arrays, 
and  use  of  adaptive  cancellers  prior  to  the  tracker  to  reject  interferences.  Sec¬ 
ondly,  it  considered  the  implementation  and  operation  of  the  tracker  in  a  sonar 
environment,  addressing  an  improved  interpolator,  multi-target  tracking,  beam- 
to-beam  handover,  and  combined  broadband  and  narrowband  tracking.  Finally,  it 
utilized  sea  tape  data  to  evaluate  tracker  performance  of  the  broadband  and  nar¬ 
rowband  trackers  with  actual  targets. 

This  report  is  organized  into  five  sections,  including  this  introduction.  Sec¬ 
tion  2. 0  discusses  the  system  aspects  of  applying  the  adaptive  tracker  in  an  actual 
sonar  environment  and  considers  the  efficient  implementation  of  the  algorithm. 

An  overall  systems  block  diagram  of  a  potential  adaptive  tracking  system  is  de¬ 
veloped.  Section  3. 0  summarizes  the  results  of  the  performance  analyses  and 
simulations  done  during  the  third  phase.  These  are  presented  without  detailed  de¬ 
scriptions  and  derivations,  which  are  included  as  appendices.  Section  4. 0  presents 
the  results  of  simulations  with  sea  tapes  of  broadband  targets.  Section  5.0  provides 
an  overview  of  the  three  phases  of  the  Adaptive  Tracking  System  Study  and  presents 
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conclusions  resulting  from  the  results  of  all  three  phase.  In  addition  to  these  five 
sections*  there  is  a  classified  supplement  to  this  report  which  deals  with  simula¬ 
tion  of  the  tracker  with  classified  sea  tapes. 
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SECTION  2.0 

SYSTEMS  CONSIDERATION  FOR  ADAPTIVE 
TRACKER  IMPLEMENTATION 


2.0  SYSTEMS  CONSIDERATIONS  FOR  ADAPTIVE  TRACKER  IMPLEMENTATION 

Most  of  the  analyses  and  simulations  of  the  adaptive  tracker  to  date  have 
attempted  to  isolate  the  sensitivities  of  the  tracker  to  various  input  parameters 
by  considering  relatively  simple  input  signals,  such  as  pure  broadband,  single 
sinusoid,  broadband  plus  a  single  sinusoid,  etc.  In  processing  sea  tape  data, 
however,  the  inputs  are  not  controlled,  and  the  tracker  must  accommodate  the 
range  of  signal  and  noise  environments  that  would  be  encountered  by  an  opera¬ 
tional  tracker.  These  include  combined  broadband  and  narrowband  targets, 
broadband  and  narrowband  interferences,  multiple  targets,  multipath  arrivals 
from  a  given  target,  and  beam-to-beam  handover  of  moving  targets.  These  issues 
are  simplified  significantly  by  first  considering  the  function  of  the  tracker  and  the 
way  in  which  it  would  be  used  operationally. 

Inherent  in  the  function  of  a  tracker  is  the  assumption  that  a  source  of 
interest  has  been  detected,  and  that  the  tracker  is  assigned,  by  an  operator,  to 
track  that  source.  For  the  broadband  tracker,  this  implies  that  the  nominal 
source  bearing  is  known,  so  that  the  split  beams  can  be  steered  such  that  the 
peak  of  the  weight  vector  appears  in  a  relatively  short  adaptive  filter.  The  actual 
detection  and  determination  of  the  approximate  bearing  will  have  been  made  by  an 
available  search  sonar.  Assignment  of  a  narrowband  tracker  requires  not  only 
detection  and  nominal  bearing  of  a  source,  but  knowledge  that  the  source  radiates 
a  narrowband  component.  This  knowledge  implies  detection  of  the  narrowband 
component  and  the  availability  of  its  approximate  frequency.  The  narrowband 
tracker  is  then  assigned  to  given  narrowband  components  by  providing  the  tracker 
with  this  coarse  frequency  estimate,  as  well  as  the  bearing  required  to  steer  the 
split  beams.  If  multiple  lines  or  multiple  targets  are  present,  the  operator 
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selects  which  of  these  is  appropriate  for  tracking,  and  assigns  the  tracker 
according  to  his  judgment.  Multiple  target  and  line  considerations  are  of  interest 
only  in  that  they  may  interfere  with  the  ability  of  the  tracker  to  maintain  track  on 
the  source  to  which  it  was  assigned.  There  is  no  need  for  the  tracker  to  determine 
the  presence  of  multiple  targets  or  to  choose  among  those  present. 

With  this  bit  of  philosophy  in  mind,  issues  of  interest  in  the  practical  imple¬ 
mentation  of  the  tracker  are  discussed  below. 

2. 1  Combined  Implementation  of  a  Broadband  and  Narrowband  Tracker 

One  of  the  attractive  features  of  the  adaptive  tracker  is  its  ability  to  track, 
via  the  weight  vector,  narrowband  and  broadband  signal  components  simultaneously. 
To  extract  the  bearing  related  to  a  broadband  component,  the  peak  of  the  time 
domain  weight  vector  is  located,  and  its  time  coordinate  converted  to  a  bearing. 

For  the  narrowband  signal,  the  weight  vector  is  Fast  Fourier  Transformed 
(FFT'd)  and  the  phase  of  the  bin  with  the  largest  magnitude  is  converted  to  a 
bearing.  This  assumes  that  the  FFT  has  sufficient  resolution  that  the  narrowband 
c  omponent  is  approximately  bln  centered. 

Two  avenues  are  possible  to  provide  increased  resolution  in  time  delay, 
hence,  increased  bearing  accuracy  for  broadband  targets;  use  of  a  high  sampling 
rate,  many  times  the  Nyqulst  rate,  or  interpolation  between  discrete  taps  at  a 
sample  rate  only  slightly  above  the  Nyquist  rate  (some  oversampling  is  required 
to  produce  correlated  samples  so  that  the  interpolator  will  work).  The  high 
computational  cost  of  the  LMS  algorithm  makes  use  of  a  high  sample  rate 
impractical.  This  combined  with  the  computational  efficiency  and  performance 
of  the  quadratic  interpolator  makes  the  choice  of  a  lower  sample  rate  obvious. 
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The  ability  to  operate  near  the  Nyquist  rate  may,  in  fact,  be  an  advantage  of  the 
adaptive  tracker  over  other  methods. 

In  the  narrowband  tracker,  frequency  resolution  serves  two  functions.  It 
provides  sufficient  gain  against  the  broadband  background  noise  to  allow  extraction 
of  a  high  quality  estimate  of  both  the  signal  phase  and  frequency  as  needed  for 
determination  of  bearing.  Analogous  to  the  broadband  case,  this  could  be  pro¬ 
vided  by  high  resolution  in  the  frequency  domain  weight  vector,  or  by  interpolating 
between  somewhat  coarser  frequency  bins.  In  simulations  to  date,  the  resolution 
has  been  provided  by  using  narrow  frequency  cells  in  the  frequency  domain  weight 
vector.  For  these  purely  narrowband  simulations,  the  total  bandwidth  has  been 
relatively  narrow,  so  this  approach  has  not  led  to  large  filter  lengths.  In  a 
system  for  combined  broadband  and  narrowband  use,  this  will  not  be  die  case, 
however.  Consider  a  system  with  total  bandwidth,  Bt,  and  sample  rate  kBt. 

For  the  sonar  to  haw  a  beamwldtb  of  Bg.  tbs  array  aperture  must  be  nominally 


where  X.  is  the  wavelength  at  the  design  frequency.  Let 


so  that  the  design  frequency  coincides  with  the  high  end  of  the  band,  so 

L=i^f-  (6) 

®BBt 
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The  maximum  propagation  time  between  array  phase  centers  (corresponding  to 
a  target  that  is  endfire  relative  to  those  phase  centers)  is  then 


'max 


L 

2c 


25 

9BBt 


(7) 


In  practice  the  actual  maximum  will  be  much  less  because  a  single  beam  will  not 
track  over  the  entire  range  of  broadside  to  endfire.  The  number  of  taps  must 
provide  ±  Tmax  seconds, 

2t  =  MTs 
max 


or 


M  =  2  Tmax  50/9BBt  =  50k 


i/kB. 


0 


B 


(8) 


If  k  is  on  the  order  of  3  to  5,  then  in  most  sonars,  less  than  50  taps  would  be 
required.  On  the  other  hand,  if  a  frequency  resolution  of  Fr  is  required, 
then 


so 

kB 

M  =^rL  (10) 

It  is  not  at  all  uncommon  for  B,/F  to  be  1000  or  more  in  narrowband  detection 

I  r 

systems,  so  clearly  if  the  filter  size  is  used  to  provide  resolution  then  the  fre¬ 
quency  domain  estimator  is  the  primary  determinant  of  M,  and  the  time  extent  of 
the  filter,  MTs,  will  be  much  larger  than  is  required  for  the  broadband  tracker. 
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The  primary  computational  load  of  the  adaptive  tracker  is  the  LMS  algorithm 
itself,  particularly  if  it  is  implemented  in  the  time  domain.  The  number  of 
computations  in  the  LMS  algorithm  is  proportional  to  the  length  of  the  adaptive 
filter,  M,  so  that  the  large  filter  sizes  needed  to  provide  frequency  resolution 
over  the  entire  input  band  will  drive  the  computational  cost  far  beyond  that 
required  by  the  broadband  tracker.  This  can  be  mitigated  somewhat  by  a 
frequency  domain  implementation  of  the  LMS  algorithm,  but  as  discussed  in  the 
following  section,  this  is  not  always  practical. 

The  size  of  the  adaptive  filter  not  only  impacts  the  computational  cost  per 
iteration  of  the  tracker.  The  algorithm  noise  is  directly  proportional  to  the  filter 
length,  so  that  providing  the  longer  adaptive  filters  necessary  for  narrowband 
tracking  will  reduce  the  performance  of  the  tracker  unless  this  noise  is  reduced 
by  decreasing  H-*  However,  decreasing  h- ,  which  Increases  the  time  constant  of 
the  filter,  degrades  the  performance  of  the  tracker  in  the  presence  of  target 
dynamics  or  changing  spectra.  Further,  use  of  small  values  of  p.  may  require 
higher  precision  in  the  implementation  of  the  device. 

Given  the  effects  of  the  filter  length,  M,  on  the  computational  cost  and 
performance  of  the  adaptive  tracker,  it  makes  sense  to  maintain  as  coarse  a 
resolution  as  possible  in  the  computation  of  the  frequency  domain  weight  vector. 
Then,  using  the  approximate  value  of  the  frequency  of  the  narrowband  component 
provided  in  the  assignment  of  the  tracker,  some  method  of  interpolation  could 
hopefully  be  used  in  that  vicinity  to  provide  higher  resolution  over  a  limited  band. 
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In  Appendix  A  of  [2],  the  mean  and  variance  of  the  frequency  domain  weight 
vector  were  derived  for  a  sinudoid  embedded  in  white  noise.  The  steady  state 
mean  weight  in  the  k^1  FFT  bin  is 


E  [  Wk(«)  J  = 


|G(k,aio)  | 2  crs2 


|G(k.u0)l 


2  2  2 

Ta  +  M  <r„ 


-iu>  —  sin  6 
e  J  o  c 


(11) 


where 


|G(k.u»0)|2 


sin2 

ft  k-  «oVl 

sin2 

i(  k-  %T.) 

(12) 


u>o  =  sinusoid  radian  frequency 

d  =  distance  between  array  phase  centers 

c  =  speed  of  sound 

0  =  bearing  of  target 

2 

<rg  =  signal  power 
<rQ2  ■  noise  power 


From  (11),  it  can  be  seen  that  once  it  has  reached  steady  state,  the  adaptive 

filter  produces  a  d.  c.  component  in  the  frequency  domain  weight  vector  for  a 

static  sinusoidal  input.  Further,  each  bin  responding  to  the  signal  has  the  correct 

phase,  *  =  (w  d/c)  sine.  This  suggests  that  additional  FFT  resolution  can  be 
o 

provided  by  processing  each  frequency  domain  weight  sequence,  Wjjm),  with  a 
d.c.  centered  FFT  bin,  that  is,  compute 
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Note  that  this  is  a  coherent  sum  in  which  real  and  imaginary  parts  are  summed 
separately  so  that  the  phase  information  is  maintained.  Due  to  the  recursive 
nature  of  the  LMS  algorithm,  the  fluctuating  part  of  the  sequence,  W^m),  Is  some¬ 
what  correlated  from  sample  to  sample,  so  the  sequence  is  not  spectrally  white. 
Consequently,  the  combination  of  an  M  point  adaptive  tracker  with  the  post 
processing  of  (13)  is  not  equivalent  to  an  MN  point  adaptive  tracker  in  terms  of 
resolution. 

This  post-processing  can  be  used  to  substantially  reduce  the  computational 
load  of  the  narrowband  tracker  as  shown  in  Figure  3.  The  adaptive  filter 
size,  M,  is  selected  to  give  reasonable  computational  requirements,  providing 
much  coarser  resolution  in  frequency  for  the  frequency  domain  weight  vector. 

Given  assignment  of  the  narrowband  tracker  to  a  given  line  component,  the  post¬ 
processing  of  (13)  is  implemented  in  the  vicinity  of  the  line,  probably  several 
FFT  bins  to  either  side  of  the  frequency  domain  weight  known  to  contain  the  signal. 
Once  sufficient  post-processing  has  been  performed,  the  W^(f)  will  be  narrow 
bandwidth,  high  SNR  FFT  bins  with  the  frequency  spacing  of  the  original  M  point 
FFT.  From  these  several  bins,  the  Wk(f)'s,  the  one  with  the  largest  magnitude 
and  the  two  adjacent  bins  are  selected.  The  phase  is  extracted  from  the  largest 
Wfe,  and  a  quadratic  interpolation  on  the  three  selected  is  used  to  estimate  the 
line  frequency,  From  Appendix  E,  where  the  quadratic  interpolator  was  used 
to  determine  delay, 
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Figure  3.  Post  Processing  Vernier  for  Narrowband  Adaptive  Tracker 
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where  u>0(£)  is  die  frequency  estimate  on  the  l  post  processor  output,  WM(f) 

is  the  processed  weight  with  the  largest  magnitude,  T  is  the  adaptive  filter  sample 

interval,  and  M  the  adaptive  filter  length.  The  bearing  estimate  is  then  given  by 


A 

where  9  is  the  phase  of  the  largest  FFT  bin. 

2.2  Time  Domain  Versus  Frequency  Domain  Implementation 

During  the  first  phase  of  the  tracker  study  [l  1 ,  a  time  and  a  frequency 
domain  algorithm  for  the  LMS  adaptive  filter  were  developed.  In  tne  frequency 
domain  version  of  an  M  tap  time  domain  filter,  the  time  series  is  passed  through 
a  M  point  FFT  and  a  single,  complex  tap  adaptive  filter  implemented  in  each 
frequency  bin.  If  the  time  domain  filter  was  real,  only  M/2  single  tap  filters 
need  be  implemented  in  the  frequency  domain  since  the  FFT  bins  are  conjugate 
symmetric.  The  time  domain  weight  vector  is  obtained  by  inverse  transforming 
the  complex  weights  across  the  M  frequency  bins. 

The  frequency  domain  algorithm  is  analytically  attractive  in  that  it  allows 
computation  of  both  the  mean  and  variance  of  the  frequency  domain  weights  under 
certain  assumptions  [1  ].  In  many  practical  cases,  this  can  be  used  to  determine 
the  statistics  of  the  time  domain  weight  vector,  which  is  not  possible  with  signal 
present  by  direct  time  domain  methods.  Even  more  significantly,  however, 
the  frequency  domain  algorithm  provides  substantial  reduction  in  computational 
cost  relative  to  the  time  domain  algorithm.  At  each  iteration,  the  frequency 
domain  algorithm  requires  one  complex  multiply  for  each  real  multiply  in  the 
time  domain,  equivalent  to  four  real  multiplies.  However,  since  only  M/2  bins 
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are  computed,  the  frequency  domain  approach  requires  nominally  twice  as  many 
multiplies  per  iteration  as  the  time  domain  algorithm.  If  the  FFT  is  computed 
without  redundancy  and  with  no  gaps  in  the  input  time  series ,  the  frequency 
domain  algorithm  iterates  only  1/M  times  as  often  as  the  time  domain  algorithm. 
Therefore,  the  M  iterations,  the  frequency  domain  algorithm  requires  approxi¬ 
mately  2/M  times  the  number  of  multiplications  of  the  time  domain  approach. 

To  this,  the  computations  of  the  two  forward  FFTs  and  the  inverse  FFT  must  be 
added.  In  fact,  if  both  time  and  frequency  domain  weights  are  needed,  the 
inverse  FFT  in  the  frequency  domain  implementation  requires  the  same  number 
of  multiplies  as  the  forward  FFT  needed  to  obtain  the  frequency  domain  weights 
from  the  time  domain.  Therefore,  letting  NT(M)  be  the  number  of  real  multiplies 
for  an  M  point  time  domain  algorithm,  the  number  of  real  multiplies  for  the 
frequency  domain  implementation,  Np,  is 

Nf  »  M  log2^  +  (Nt(M»  (16) 

This  produces  sizeable  reductions  for  even  modest  M. 

Unfortunately,  the  frequency  domain  algorithm  is  not  an  exact  realization 
of  the  time  domain  LMS  algorithm.  It  is  subject  to  windowing  effects  due  to 
the  block  processing  in  the  FFT.  For  a  broadband  signal,  this  effect  can 
only  be  regarded  as  insignificant  if  [  2  ] 

M>>A  (17) 

where  A  is  the  delay  in  the  signal  between  adaptive  filter  inputs  in  sample 
intervals.  In  a  broadband  tracker  implemented  in  the  time  domain,  it  is  desir¬ 
able  to  keep  M  small  and  place  the  peak  near  the  center  of  the  adaptive  filter 


when  the  target  is  on  the  split  beam  MRA,  so  that  A/M=l/2.  If  this  tracker  were 
to  be  implemented  in  the  frequency  domain,  M  would  have  to  be  increased  sub¬ 
stantially  to  reduce  windowing  effects,  offsetting  some  of  the  computational 
advantages  of  the  frequency  domain  algorithm,  and  increasing  algorithm  noise. 

On  the  other  hand,  for  a  combined  narrowband  and  broadband  system  that  uses 
large  M  to  achieve  frequency  resolution,  by  the  results  of  the  preceding  section, 
the  total  time  delay  needed  for  the  broadband  signal  is  a  small  fraction  of  the 
filter  length.  By  placing  the  peak  near  the  front  of  the  adaptive  filter  for  all 
bearing  angles,  it  can  then  be  assured  that  ( 17)  is  true,  so  window  effects  will  be 
negligible. 

It  has  been  shown  that  the  variance  of  the  time  domain  weight  vector  obtained 
via  the  two  implementations  is  the  same  for  broadband  inputs  at  low  signal-to- 
noise  ratio  whenever  M  >  >  2.  Simulations  appear  to  support  this  at  higher 
signal  to  noise  ratios  as  well.  However,  it  is  possible,  given  a  particular 
desired  time  constant,  for  the  filter  to  be  stable  in  the  time  domain  but  unstable 
in  the  frequency  domain.  The  time  constant,  in  iterations  of  the  time  domain 
filter  is 

Ttl  =  2HPiQ  (18) 

where  Pin  is  the  input  power.  Due  to  the  gain  of  the  F FT,  the  time  constant  of 
the  frequency  domain  filter  is 

Tfl  =  2RMPln  (l9) 
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where  it  is  assumed  that  the  input  power  is  primarily  broadband.  Since  the 
frequency  domain  implementation  iterates  only  l/M  times  as  often,  the  two 
filters  have  the  same  time  constant  in  seconds.  However,  the  time  constant  for 
the  frequency  domain  implementation  is  much  smaller,  i.  e. ,  closer  to  the 
unstable  region.  If,  for  example,  the  system  required  t  =  256  and  M  =  1024, 
then  =  1/4  for  equal  response  in  seconds,  and  the  frequency  domain 
implementation  is  unstable. 

The  choice  of  frequency  or  time  domain  implementation  must  therefore  be 
made  on  a  case  by  case  basis,  taking  into  account  the  above  considerations. 

For  broadband  only  trackers,  the  increased  filter  size  necessary  to  use  a 
frequency  domain  implementation  will  offset  some  of  Its  computational  advantages. 
For  purely  narrowband  applications ,  the  frequency  domain  approach  should  be 
used  if  the  combination  of  M  and  the  desired  time  constant  do  not  lead  to  unstable 
or  near  unstable  operation.  Implementation  of  the  tracker  for  combined  broad¬ 
band  and  narrowband  targets  must  consider  both  factors. 

2.3  Tracking  of  Multiple  Targets 

Typically,  the  sonar  environment  will  consist  of  a  number  of  directional 
sources,  some  of  which  will  be  targets  of  interest,  others  interferences.  Often, 
the  nature  of  some  contacts  will  be  unknown,  and  targets  of  interest  will  be 
differentiated  from  interferences  only  after  tracking  for  some  period  of  time. 

The  adaptive  tracker  structure  must  therefore  have  the  capability  of  tracking  in 
this  multi-contact  environment. 

At  the  outset,  multiple  contact  tracking  is  simplified  by  the  spatial 
response  of  the  split  beams.  Given  two  contacts  of  similar  strength  that  are 
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separated  in  bearing  by  more  than  a  beamwidth,  the  output  of  a  beam  steered  at 
one  can  essentially  be  regarded  as  containing  the  single  contact,  since  the  other 
is  in  the  sidelobes  and  severely  attenuated.  Therefore,  a  tracker  structure 
operating  on  the  largest  peak  or  largest  line  in  the  weight  vector  will  almost 
always  track  the  target  of  interest .  Multiple  tracking  will  then  consist  of 
assigning  a  second  split- beam  pair  and  adaptive  tracker  to  the  other  contact. 
Hence,  for  contacts  of  approximately  the  same  radiated  noise  level  separated  by 
more  than  a  beamwidth,  the  multiple  target  tracking  problem  reduces  to  multiple 
single-target  trackers. 

The  problem  is  complicated  somewhat  when  one  contact  is  much  stronger 

than  the  other,  and  it  is  desired  to  track  the  weaker  of  the  two.  Even  though  the 

stronger  contact  is  in  the  sidelobes  of  a  beam  steered  at  the  weaker  one,  it  may 

be  sufficiently  strong  to  produce  a  peak  of  nearly  equal  or  greater  amplitude  in 

the  weight  vector.  If  the  adaptive  tracker  merely  picks  the  largest  peak,  then  the 

tracker  will  operate  on  the  target  that  is  in  the  sidelobes  whenever  its  amplitude 

exceeds  that  of  the  target  of  interest.  This  problem  can  be  eliminated  by 

restricting  the  choice  of  the  largest  weight  to  that  range  of  delays  that  would 

occur  for  targets  in  the  mainlobe,  say  ±  T _ ,  where 

max 

T  =  —  sin-*-  (20) 

max  c  2  y  ' 

with  0Q  the  beamwidth  of  the  split  beams.  This  is  most  efficiently  done  from  a 
computational  point  of  view  by  restrldtlng  the  length  of  the  adaptive  filter  to  2T 

ulaA 

and  providing  bulk  delay  to  place  a  target  on  the  beam  maximum  response  axis 
(MRA)  at  the  center  of  the  weight  vector,  as  shown  in  Figure  4. 
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In  implementations  for  combined  broadband  and  narrowband  tracking  (see 
Section  2.2),  the  adaptive  filter  length  may  be  required  to  provide  frequency 
resolution.  In  that  case,  the  restriction  to  ±T  is  made  by  limiting  the 

UlSX 

taps  used  in  the  broadband  estimator  rather  than  shortening  the  time  domain 
weight  vector.  Note,  however,  that  the  presence  of  the  strong  contact  may 
still  bias  the  location  of  the  peak,  or  even  mask  its  location.  This  situation  is 
best  handled  by  use  of  a  sidelobe  canceller  prior  to  the  tracker.  The  stronger 
contact  is  then  treated  as  an  interference,  and  handled  as  in  Section  3.5 

If  the  target  is  narrowband,  the  estimate  of  bearing  is  based  upon  the  phase 
of  a  line  component  radiated  by  the  contact.  Restriction  of  the  length  of  the 
adaptive  filter  will  not  change  the  effects  of  a  strong  contact  in  the  sidelobes. 
In  assigning  the  narrowband  tracker,  though,  a  particular  frequency  will  have 
been  selected  for  tracking  so  that  the  narrowband  estimate  can  be  restricted  to 
the  vicinity  of  that  line  in  the  FFT  of  the  weight  vector.  Unless  a  line  radiated 
by  the  sidelobe  contact  is  close  to  the  chosen  line  from  the  contact  of  interest, 
this  will  allow  the  tracker  to  operate  on  the  weaker  contact.  This  should  be 
taken  into  account  in  the  choice  of  a  line  for  assignment  to  the  tracker.  A  line 
from  a  sidelobe  contact  that  is  too  close  to  a  line  of  interest  may  be  treated  as 
an  interference  and  rejected  by  sidelobe  cancellation  (see  Section  3.5).  A 
strong  broadband  contact  in  the  sidelobes  may  bias  the  narrowband  estimate, 
as  discussed  in  [2],  but  the  sidelobe  cancellation  methods  can  also  be  used  to 
eliminate  the  effects. 

The  most  severe  case  of  multiple  target  tracking  occurs  when  multiple 
targets  appear  within  the  mainlobe  of  a  single  beam.  This  will  occur  most  often 
when  two  targets  originally  in  separate  beams  cross  in  bearing,  and  therefore 


appear  in  the  same  beam  for  some  period  of  time.  The  problem  is  then  twofold; 
maintain  contact  with  both  targets  while  they  are  in  the  same  beam,  and  dis¬ 
criminate  between  the  two  targets  so  that  the  individual  tracks  may  be  maintained. 
The  effects  of  such  an  encounter  depends  upon  whether  broadband  or  narrowband 
tracking  is  being  used  on  each  of  the  two  contacts. 

a.  Both  Narrowband  Targets:  If  the  lines  used  to  track  the  two  targets 
are  sufficiently  far  apart  in  frequency  with  respect  to  the  resolution 
of  the  frequency  domain  weight  vector,  then  the  two  tracks  can  be 
maintained.  However,  if  their  frequencies  are  indistinguishable 

by  the  algorithm,  they  will  appear  as  a  single  target  while  they  are 
in  the  main  lobe. 

b.  Broadband  Target  and  Narrowband  Target;  It  is  unlikely  that  the 
narrowband  target  would  be  sufficiently  strong  to  appreciably  affect 
the  peak  of  the  time  domain  weights,  so  it  will  be  possible  to  main¬ 
tain  track  on  the  broadband  target  with  little  degradation.  As  shown 
in  Appendix  D  of  [  2] ,  the  broadband  target  will  bias  the  narrowband 
estimate  and  increase  its  variance.  If  the  broadband  target  were 
very  strong,  the  narrowband  estimator  could  track  the  broadband 
target  instead  of  the  line  to  which  it  was  assigned. 

c.  Both  Broadband  Targets:  When  the  peaks  in  the  weight  vector 
corresponding  to  each  target  are  close  together,  their  locations 
become  biased  toward  the  other  target,  eventually  becoming  a 
single  peak,  indistinguishable  from  a  single  target  without  the  use 
of  prior  history  of  the  tracks. 

These  considerations  are  very  much  the  same  as  those  involved  in  the  use  of 
conventional  trackers  in  a  multi-target  environment. 


In  conventional  broadband  trackers ,  the  time  delay  estimate  is  generally 
based  upon  the  determination  of  the  location  of  the  peak  of  the  cross-correlation 
function  between  the  split  array  outputs  (or,  alternatively,  the  zero  crossing  of 
its  derivative).  The  correlation  function  is  computed  at  a  single  value  of  delay 
and  r  varied  to  maximize  the  cross-correlation  function.  This  is  done  as  shown 
in  Figure  5,  in  which  a  variable  delay  is  inserted  in  one  line  of  the  tracker 
to  control  the  point  at  which  the  derivative  of  the  cross-correlation  function  is 
computed.  The  sign  of  the  derivative  determines  the  direction  in  which  the  value 
of  the  current  t  is  incremented,  with  equilibrium  occurring  when  t  is  at  the  zero 
crossing.  Typically,  the  output  of  the  correlator,  y(n),  is  fed  back  through  some 
type  of  tracking  logic  which  controls  how  t  is  modified  for  a  given  value  of 
y(n).  In  most  systems,  the  tracking  logic  includes  smoothing,  and  a  first  order 
update  of  the  delay,  that  is 

t  (n+i)  =  r(n)  -  ®Ys(n)  (20) 

where 

T(n)  =  nth  value  of  delay 
and 

Yg(n)  =  nth  smoothed  value  of  y 

Higher  order  tracking  loops  vary  T(n)  based  on  not  only  the  value  of  y(n),  but  its 
derivatives  (actually,  higher  order  differences  in  a  discrete  system). 

The  incorporation  of  such  loops  into  the  tracker  provides  the  tracker  with 
the  ability  to  follow  dynamic  signals.  For  example,  the  first  order  loop  can 
follow  a  linearly  changing  delay  with  a  fixed  lag  error  (the  lag  being  dependent 
upon  the  gains  incorporated  in  the  loop).  However,  the  tracking  loop  serves 


another  function  whenever  two  targets  cross  in  bearing,  that  of  allowing  the 
tracks  to  be  maintained  using  the  past  dynamics  of  the  target.  Clearly,  if  one  of 
the  targets  changes  its  dynamics  during  the  time  in  which  the  zero  crossings  in 
the  derivative  are  indistinguishable,  then  loss  of  track  may  occur. 

A  similar  approach  can  be  applied  to  the  adaptive  tracker  to  handle  the  case 
when  two  broadband  targets  appear  in  one  beam.  A  relatively  narrow  window  is 
placed  about  the  current  peak  in  the  weight  vector,  and  the  peak  selection  is  only 
allowed  to  occur  within  the  window.  For  example,  if  the  current  delay  estimate 
from  the  interpolator  is  9  (m),  then  on  the  (m+l)th  iteration  the  peak  (which 
provides  the  starting  point  for  the  interpolator)  will  be  the  largest  weight  within 
±^max  °f  *(m).  This  prevents  the  interpolator  from  jumping  to  another  target 
more  than  Tmgy  from  the  current  location.  Clearly,  if  the  window  is  too  narrow 
(Tmax  400  sma^)  then  the  dynamics  of  the  source  can  carry  the  peak  out  of  the 
window  between  iterations,  and  track  will  be  lost,  so  the  window  width  must  take 
into  account  die  maximum  rate  of  change  of  delay. 

If  the  current  estimate,  T(m),  were  used  to  control  the  next  location  of  the 
window,  then  whenever  the  fluctuations  in  T(m)  were  large  in  comparison  to 
Tmax'  tben  a  sin8le  fluctuation  could  move  the  window  away  from  the  peak,  and 
again  track  would  be  lost.  This  can  be  prevented  by  passing  the  estimate  T(m) 
through  tracking  logic  prior  to  moving  the  window,  just  as  is  done  in  conventional 
trackers.  The  width  of  the  window  then  must  consider  the  signal  dynamics  and 
the  fluctuations  in  the  output  of  the  tracking  logic.  The  simplest  tracking  logic 
would  simply  smooth  the  sequence  T(m),  but  in  order  to  maintain  track  using 


target  dynamics,  higher  order  differences  of  the  smoothed  delay,  T.(m),  are 

s 

used  to  update  the  current  window  location,  as  shown  in  Figure  6.  A  first 
order  loop  on  the  window  would  compute  the  new  window  location  as 

Tw(m)  =  Tg(m)  -  otTw<m-l) 

for  example.  Now,  when  two  windows  overlap,  only  a  single  estimate  would  be 
supplied  to  the  loops  controlling  the  windows  analogous  to  the  conventional 
tracking  loops.  If  die  targets  do  not  execute  maneuvers  during  the  period  when 
the  windows  overlap,  the  loop  will  allow  maintenance  of  the  two  tracks. 

Use  of  the  loops  to  control  the  location  of  windows  operating  on  the  weight 
vectors  instead  of  varying  a  bulk  delay  prior  to  the  input  of  the  adaptive  filter 
has  the  advantage  of  decoupling  the  tracking  loop  from  the  dynamic  behavior 
of  adaptive  filter.  If  the  loop  controlled  a  delay  prior  to  the  tracker,  as  shown 
in  Figure  6,  then  the  dynamic  behavior  of  the  combination  of  the  tracking 
loop  and  adaptive  filter  (which  exhibits  its  own  first  order  behavior  in  the  weight 
updates)  would  be  extremely  difficult  to  quantify  given  the  nonlinear  properties 
of  the  adaptive  filter.  By  using  the  windowing  scheme,  the  dynamics  and  statistics 
of  the  adaptive  tracker  estimate  can  be  characterized  using  the  methods  developed 
in  [1]  and  [2] ,  then  the  response  of  the  loop  to  these  input  statistics  analyzed 
by  linear  system  theory.  This  of  course  assumes  that  the  window  is  sufficiently 
wide  that  it  does  not  affect  the  statistical  behavior  of  the  interpolator.  It  has  been 
shown  in  (21  that  the  adaptive  tracker  can  follow  a  linearly  changing  delay  with 
a  fixed  lag.  Therefore,  a  first  ojrder  loop  would  simply  introduce  an  additional 
lag  in  the  window  location.  Taking  this  into  account  in  the  width  of  the  window 
would  allow  the  estimator  to  track  a  ramp  input  without  any  additional  lag,  since 


27 


Figure  6.  Delay  Windowing  with  a  Tracking  Loop 


the  window  lag  does  not  affect  the  estimate  (unless  it  becomes  too  great,  so  that 
track  is  lost).  Use  of  a  second  order  loop  would  allow  tracking  of  a  ramp  input 
without  accounting  for  any  lag  in  the  window  width,  although  some  accounting  for 
the  increased  fluctuations  in  the  current  window  center  would  be  necessary.  It 
must  be  noted  that  because  the  behavior  of  the  adapter  filter  and  the  window 
tracking  loop  are  decoupled,  the  order  of  the  combination  is  limited  by  that  of  the 
adaptive  tracker.  That  is,  use  of  a  higher  order  loop  in  the  windowing  does  not 
give  higher  order  behavior  for  the  tracker  since  the  adaptive  filter  is  first  order. 

2.4  Beam-to-Beam  Handover  of  Dynamic  Targets 

Whenever  split  array  beamformers  are  used  to  provide  the  input  to  the 
adaptive  tracker,  the  potential  exists  for  moving  targets  to  move  out  of  the  main 
lobe  of  the  split  beam,  resulting  in  loss  of  track.  Two  methods  of  maintaining 
track  are  steering  of  the  split  beams  to  maintain  contact ,  or  handover  of  the 
tracking  function  to  an  adjacent  fixed  beam.  Steering  of  the  split  beams 
requires  a  beamformer  structure  that  is  dedicated  to  the  tracker  and  capable 
of  relatively  high  resolution  steering.  It  also  requires  incorporation  of 
a  tracking  loop  to  control  the  steerable  beam.  Beam-to-beam  handover  applies 
when  the  tracker  is  used  as  an  add-on  to  a  preformed  beam  sonar,  and  minimizes 
the  cost  of  the  tracker  itself.  This  is  the  configuration  that  is  considered  here. 

Suppose  a  dynamic  target  is  passing  through  a  beam  to  which  the  adaptive 


tracker  is  assigned.  The  goal  of  the  beam-to-beam  handover  is  to  reassign 
the  tracker  to  the  next  beam  before  it  is  sufficiently  attenuated  to  lose 
track.  Since  the  array  pattern  is  known,  the  bearing  angle  corresponding 
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to,  say,  a  X  dB  attenuation  can  be  determined  as  a  potential  handover  point. 

When  the  bearing  estimate  reaches  this  point,  then  handover  is  initiated.  Several 
factors  must  be  taken  into  account  in  implementation  of  the  handover,  however; 
a.  The  bearing  estimate  from  the  tracker  is  random,  and  as  such, 
subject  to  fluctuations.  When  the  target  is  in  the  vicinity  of  the 
handover  point,  the  estimate  may  fluctuate  back  and  forth  about 
the  desired  point.  To  avoid  rapid  handover  of  the  target  back  and 
forth  during  this  time,  the  handover  algorithm  must  include  some 
hysteresis.  The  simplest  means  of  doing  this  is  to  designate  a 
different  value  for  handover  when  the  target  passes  to  or  from  a 
beam.  For  example,  let  4>i  and  4>.+1  be  the  location  of  two  adjacent 
beams,  and  let  the  desired  handover  point  be 


hd  - 


*i+Vi 


(22) 


Then  when  a  target  passes  from  beam  i  to  beam  i+1,  the  actual 
handover  is  performed  when  the  bearing  estimate,  4>,  reaches 

hd  (i,  i+1)  =  hd  +  6d  (23) 

while,  when  it  passes  from  beam  i+1  to  beam  i,  it  occurs  when  <t> 
reaches 

hd(i+l,i)  «  hd  -  6d  (24) 
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There  is  therefore  a  "dead  zone"  of  26^  degrees  over  which  the 
estimate  can  fluctuate  after  beam  handover  that  will  not  initiate 
another  handover  back  to  the  original  beam.  If  6^  is  sufficiently 
large  with  respect  to  the  fluctuations  in  4>,  the  problem  of  jittering 
at  handover  will  be  eliminated. 

b.  Since  the  target  is  dynamic,  the  estimate  <b  will  lag  the  actual 
bearing  of  the  target,  4>,  so  that  handover  will  actually  occur  later 
than  desired  if  <t>  is  used  to  initiate  the  process.  This  is  compounded 
by  the  fact  that  incorporation  of  the  "dead  zone"  in  (a)  makes  hand¬ 
over  occur  later  than  the  desired  point,  h^,  anyway.  In  many  sonars, 
the  beam  spacing  will  be  sufficiently  fine  (i.  e. ,  beam  crossover 
points  not  highly  attenuated)  compared  to  lag  that  this  effect  will  not 
be  a  problem.  If  the  lag  is  a  problem,  a  possible  solution  is  to 
estimate  the  lag  from  the  bearing  rate  and  add  it  to  4>  prior  to 
comparison  to  the  beam  handover  point. 

c .  The  beam  handover  process  should  include  the  handover  of  bearing 
information  from  the  original  beam  to  the  new  beam.  In  the 
adaptive  tracker,  this  information  is  in  the  adaptive  filter  weight 
vector.  The  actual  shape  of  the  weight  vector  should  not  change 
drastically,  since  the  spectrum  of  the  source  remains  the  same, 
with  any  changes  due  to  the  difference  in  the  frequency  response 

of  the  beams.  The  primary  change  will  be  in  the  form  of  a  delay  in 
the  weights  corresponding  to  the  change  in  the  delay  between  phase 
centers  for  the  two  beams.  Because  the  phase  centers  are  known 
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for  each  of  the  fixed  beams,  hence,  their  separation,  dj,  the  time 
delay  for  the  new  beam  can  be  estimated  from  the  bearing  estimate. 
Let  the  index  k  denote  the  beam  from  which  the  target  is  being 
handed  and  k+1  denote  the  new  beam,  and 
=  bearing  of  beam  j 

dk+l  =  effective  distance  between  phase  centers  for  beam  j 
xk+1  =  estimate  of  delay  between  phase  centers  for  beam  j 
<p  =  current  bearing  estimate 
Then 

A  ^k+l 

Tk+1  =  —  sin  <♦  -  *k+1) 

where  c  is  die  speed  of  sound.  The  weight  vector  for  the  new  beam 
is  therefore  approximately  the  weight  vector  for  the  original  beam 
shifted  by  ($k+1  -  t),  which  would  be  quantized  to  the  sample  rate 
for  implementational  simplicity. 


Two  alternatives  exist  for  inserting  this  shift,  actually  time  shifting 
the  time  domain  weight  vector  or  by  incorporating  a  phase  shift 
corresponding  to  this  delay  in  each  FFT  bin  of  the  frequency  domain 
weight.  Implemented  in  the  time  domain,  the  initial  weight  vector 
for  the  new  beam  would  be 


Wk+i(i)  =  Wk(i 


A  A 

X  -  T 

A  * 

x  _x 

k+1  k 

,  ISi  - 

k+1  k 

T,  J 

l  Ts 

<  M 


(26) 


where 

a  *h 

Wk+1  =  i  initial  weight  for  new  beam 

Wk(i)  =  ltd  weight  from  original  beam 
tXj  a  largest  integer  £  X 
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This  leaves  M  - 


weights  undefined.  Possible  alternatives 


?k+l 

[  Ts  . 

for  these  weights  are  an  end  around  shift  of  the  weight  vector  or 
arbitrary  assignment  of  these  weights  to  zero.  If  the  length  of  the 
adaptive  filter  is  judiciously  chosen,  the  shift  will  be  valid  in  the 
vicinity  of  the  peak  of  the  time  domain  weight  vector,  hence 
acceptable  in  terms  of  the  broadband  estimator.  However,  either 
end  around  shift  or  arbitrary  zero  assignment  of  the  weights  near 
the  end  of  the  filter  will  reduce  the  coherence  of  any  narrowband 
components.  The  end  around  shift  will  minimize  this  degradation. 


.th 


If  Fk(i)  is  the  i  frequency  domain  weight  of  the  original  beam 

and  Fk+1(i)  is  the  1th  initial  frequency  domain  weight  for  the  newest 

beam,  then  the  phase  shift  method  gives 

.gjLi  |Vi  :!ki 

'Ml  T  J 


W‘>  ’  V* e 


(27) 


This  produces  exactly  the  end  around  time  shift  described  above. 
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3.0  SUMMARY  OF  PERFORMANCE  ANALYSIS 

This  section  summarizes  the  results  of  the  performance  analyses  performed 
during  the  third  phase  of  the  Adaptive  Tracking  System  Study.  Section  2. 1  con¬ 
siders  the  effects  of  multipath  on  the  adaptive  tracker,  particularly  as  they  re¬ 
late  the  multipath  effects  in  conventional  trackers.  Section  2. 2  discusses  the  ef¬ 
fects  of  arrays  on  tracking  of  both  static  and  dynamic  targets.  This  extends  results 
of  the  first  two  phase  of  this  study,  which  considered  each  half  of  the  split  array  as 
a  point  omnidirectional  hydrophone. 

In  Section  2. 3,  estimators  are  developed  for  tracking  of  targets  which  have  both 
broadband  and  narrowband  components.  Several  structures  are  developed  which 
require  differing  amounts  of  a  priori  knowledge  of  the  signal  and  noise  character¬ 
istics.  Section  2.4  described  a  simplified  interpolator  for  the  broadband  tracker, 
and  analyzes  its  performance  relative  to  the  sin  x/x  interpolator  used  in  Cl 3. 

Section  2. 5  analyzes  the  effects  of  using  adaptive  cancellers  on  the  split  beam 
outputs  prior  to  tracking.  Finally,  Section  2. 6  summarizes  the  results  of  simula¬ 
tions  of  the  adaptive  tracker  and  a  conventional  tracker  with  moving  targets. 

3. 1  Effects  of  Multipath  on  the  Adaptive  Tracker.  The  presense  of  multipath 
propagation  causes  multiple  delayed  and  attenuated  replicas  of  the  target  wave¬ 
form  to  appear  at  the  split  beam  output.  It  is  shown  in  Appendix  A  that  these  mul¬ 
tiple  arrivals  can  result  in  displacement  and  broadening  of  the  peak  of  the  adaptive 
tracker  weight  vector  for  broadband  targets,  and  may  produce  ambiguous  peaks. 
These  effects  will  depend  upon  the  particular  propagation  path  lengths  involved  and 
upon  the  arrival  angles  in  both  bearing  and  depression/elevation  (d/e).  The  poten¬ 
tial  effects  of  multipath  therefore  include  bias  in  the  estimate  due  to  displacement 
of  the  peak  or  selection  of  an  ambiguous  peak,  and  increased  variance  due  to  peak 
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broadening.  These  effects  are  not  readily  generalized  since  they  apply  to  each 
specific  multipath  structure.  However,  it  is  shown  that  the  broadband  adaptive 
tracker  is  affected  by  multipath  in  the  same  way  as  conventional  trackers. 

Similarly,  it  is  shown  that  multipath  produces  multiple  sinusoids  of  differing 
phase  in  the  weight  vector.  In  dynamic  situations,  it  may  be  possible  to  resolve 
these  multiple  sinusoids  by  means  of  differing  doppler  of  different  vertical  arrival 
angles.  Such  lines  produce  ambiguities  analogous  to  multiple  peaks  in  the  broad¬ 
band  tracker.  Lines  that  are  not  resolvable  by  the  narrowband  tracker  will  pro¬ 
duce  a  biased  phase  estimate  from  the  FFT  bin,  and  will  cause  increased  variance 
due  to  interference  of  the  sinusoids.  Again,  the  quantitative  effects  depend  upon 
the  specific  multipath  structure  present,  but  it  is  shown  that  the  adaptive  tracker 
is  affected  In  the  same  way  as  conventional  trackers. 

3.2  Effects  of  Arrays  on  the  Adaptive  Tracker.  Throughout  the  first  two 
phases!1 » of  this  study,  the  inputs  to  the  adaptive  tracker  have  been  assumed 
to  come  from  two  omnidirectional  hydrophones  separated  by  a  distance,  d.  It  was 
assumed  that  arrays  could  replace  the  omnidirectional  hydrophones  with  minimal 
modification  to  the  results.  This  section  considers  the  effects  of  such  arrays  on 
previous  analyses  in  both  the  static  and  dynamic  cases. 

In  the  static  case,  it  is  shown  in  Appendix  B  that  the  addition  of  arrays  pro¬ 
duces  the  equivalent  model  of  the  adaptive  tracker  shown  in  Figure  7.  The  results 
developed  for  the  omnidirectional  hydrophones  can  therefore  be  extended  by  re¬ 
placing  the  half  arrays  and  beamformers  with  linear  filters,  Hx(o>)  and  H^w),  1.  e. , 
simply  by  modifying  the  Input  spectra  in  previous  results.  Most  significantly,  it  is 
shown  that  when  the  two  split  arrays  are  Identical  with  identical  shadings  the  equiv¬ 
alent  model  is  as  in  Figure  8.  Therefore,  the  analyses  using  omnidirectional 
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hydrophones  in  Cl]  and  C2]  can  be  applied  when  arrays  are  used  with  static  tar¬ 
gets  by  modifying  the  signal  spectrum  by  the  transfer  function,  Hx(u>),  and  letting 
d  be  the  distance  between  split  array  phase  centers. 

Analysis  of  the  dynamic  case,  in  which  the  delay  associated  with  the  target 
changes  linearly  with  time,  is  complicated  by  the  fact  that  the  target  will  eventually 
leave  the  main  beam,  making  steady  state  analysis  impossible.  This  was  overcome 
by  making  the  time  delay  periodic  across  the  main  beam,  performing  a  steady 
state  analysis,  then  identifying  and  removing  terms  associated  with  the  periodicity 
in  the  results.  It  is  shown  that  the  lag  in  the  peak  of  the  weights  at  low  SNR,  cor¬ 
responding  to  the  lag  in  the  bearing  estimate  is 

Lag  m  — j  0.69  +  21nMf  (28) 

(pi<rn  -  brd)d  M 

where 

M  =  number  of  hydrophones  in  half  array 
r  =  linear  rate  of  change  of  delay 
6  3  algorithm  sample  rate 
b  =  bandwidth  of  input  processes 
H  3  feedback  coefficient  of  the  adaptive  algorithm 

<r^  3  noise  power. 

This  agrees  with  the  omnidirectional  hydrophone  case  derived  in  C2]  when  M*  1. 
Asymptotically,  the  lag  is  directly  proportional  to  InM/M,  thus  decreasing  with  the 
number  of  array  elements. 

3.3  Estimator  Structures  for  Use  With  Both  Broadband  and  Narrowband 
Components.  As  pointed  out  in  the  introduction,  one  of  the  motivations  for  use  of 
an  adaptive  tracker  is  that  it  has  the  potential  to  track  signals  using  both  broadband 
and  narrowband  energy  simultaneously.  Reference  Cl 3  developed  estimator 


structures  for  use  when  either  broadband  or  narrowband  energy  only  were  present 
In  the  radiated  noise  of  the  target.  References  Cl]  and [ 2]  then  analyzed  the  per¬ 
formance  of  these  two  estimators  in  the  presence  of  targets  that  are  broadband, 
narrowband,  or  a  combination  of  the  two.  However,  a  means  of  combining  the  two 
estimates,  or  of  performing  a  different  estimate  when  the  target  has  both  broad¬ 
band  and  narrowband  components  was  not  considered. 

In  Appendix  D,  three  estimator  structures  are  derived  for  use  when  both 
broadband  and  narrowband  components  are  present  in  the  target  spectrum.  These 
three  structures  require  differing  amounts  of  a  priori  knowledge  as  to  the  signal. 

It  may  readily  be  argued  that  the  weights  of  the  adaptive  filter  are  gausslan, 
and  an  expression  for  the  mean  weights  with  combined  narrowband  and  broadband 
target  was  derived  in  C 23.  Further,  the  results  of  C 2D  show  that  the  weight  vari¬ 
ance  is  Independent  of  the  bearing.  With  the  additional  (approximate)  assumption 
that  the  weights  are  uncorrelated,  it  is  possible  to  determine  the  maximum  likeli¬ 
hood  (ML)  estimator  for  bearing  given  the  adaptive  filter  weights.  If  the  broadband 

f 

and  narrowband  SNR's  are  also  unknown  then  they  are  jointly  estimated  in  the  ML 
structure.  The  ML  estimator  maximizes  the  functional 

J(fl)  *  Qx  (£)  +  Q2  f2  (3)  +  Q2  ($>  f2  ($)  (29) 

where 

(9)  -  |wf(max)  |  cos  (*w  -  wq|  sin  ?)  (30) 

and 

.  M-i  /  A  /N\ 

f2  (0)  *  Z  W(k)  p  (kTs  -  S  sin  9)  (31) 

k*0  \  c  / 
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and  where  Q^,  Q2,  and  Q3  are  constants  essentially  independent  of  0,  given  in  Ap¬ 
pendix  D.  Also 

|  W{(max)  |  =  magnitude  of  the  frequency  domain  weight  with  the  largest 
magnitude 

=  phase  of  the  frequency  domain  weight  with  the  largest  magnitude 
W(k)  =  time  domain  weight 

p(kTs)  =  interpolation  function  with  the  shape  of  the  broadband  mean 
weights 

M  *  number  of  taps  in  the  adaptive  filter. 

This  estimation  structure  is  shown  in  Figure  9. 

If  the  broadband  and  narrowband  signal  to  noise  ratios  are  known  a  priori, 
or  can  be  estimated  from  the  sonars,  then  the  estimator  structure  simplifies  con¬ 
siderable.  The  ML  estimate  then  maximises 


J^) 


M7NB +  VBB 


+  1  (Ybb  +  1)  fl(?)  + 


-lag-  f  ( 

VBB  +  l  2< 


where  and  ?BB  are  the  narrowband  and  broadband  signal  to  noise  ratios, 
respectively.  This  structure  is  shown  in  Figure  10. 

An  attractive  feature  of  these  structures  is  that  they  retain  elements  of  the 
estimators  developed  in  [2]  for  purely  broadband  or  narrowband  signals.  In  fact, 
the  second  estimator  reduces  the  appropriate  estimator  whenever  either  or 
?Bb  18  zero*  This  is  a  strong  justification  of  the  original  estimators. 

Appendix  D  also  presents  an  alternate  structure  which  requires  no  a  priori 
knowledge  of  the  signal  and  noise  characteristics,  but  combines  the  broadband  and 
narrowband  estimates  based  upon  their  estimated  variance.  The  combined  estimate 
is 
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MINIMIZATION  LOOP 

Figure  9.  Combined  Adaptive  Tracking  System  with  Unknown  SNRs 
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where 
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rl 


=  estimated  variance  of  the  narrowband  estimate 


estimated  variance  of  the  broadband  estimate 
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current  narrowband  estimate 
current  broadband  estimate. 


The  values  of  a\  and  <r%  are  obtained  by  moving  averages  over  and 
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If  the  estimator  variances  were  known  exactly,  and  and  were  uncor- 

relfted,  this  structure  would  produce  a  combined  estimate  with  smaller  variance 

than  that  of  or  f__.  The  fluctuations  in  the  estimates  of  of  and  <r\  will 
BB  ^NB  ^BB 

Increase  this  variance.  Further,  correlation  between  3^  and  ^Bfi  will  increase 

the  variance  of  -the  estimate.  Appendix  O  includes  an  analysis  of  the  sensitivity  of 

this  combiner  to  correlation. 

3.4  Analysis  of  a  Simplified  Interpolator.  The  broadband  adaptive  tracker, 
as  originally  described  In  C  1 3,  used  a  sin  x/x  interpolator  to  detennine  the  loca¬ 
tion  of  the  peak  In  the  weights  when  it  fell  between  discrete  taps.  This  procedure 
was  computationally  inefficient  in  that  it  utilised  all  tracker  weights  and  required 
an  iterative  procedure  to  find  tbs  peak.  Since  many  of  the  weights  are  dominated 
by  fluctuations,  use  of  all  the  weights  is  also  a  questionable  procedure  from  a  per¬ 
formance  point  of  view. 
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In  Appendix  E,  a  method  of  locating  the  peak  based  upon  a  quadratic  fit  to  the 


largest  weight  and  the  two  adjacent  weights  is  described.  This  interpolator  has 
the  computational  advantage  that  the  peak  is  computed  directly  from  the  values  of 
the  largest  weight  and  the  two  adjacent  weights  as 


wj-i  - 

wj-i  "  wj+i  - 2' 


where  Wj  is  the  largest  weight,  the  adjacent  weight  with  smaller  time  index, 

and  WJ+1  the  other  adjacent  weight.  Also,  J  is  the  time  index  of  the  largest  weight 

and  T  the  algorithm  sample  rate.  Note  that  no  Iterative  procedure  is  needed  to 
s 

find  the  peak. 

The  variance  of  this  estimate,  r,  is  then  determined  analytically  using  the 
method  used  to  determine  the  variance  in  peak  localization  in  C  l3  and  C  23,  origin¬ 
ally  described  in  C43.  These  results  are  compared  to  the  variance  of  the  estimate 
using  the  sin  x/x  interpolator  developed  in  Appendix  V  of  [  1  ].  The  results  of  the 
comparison  must  be  evaluated  numerically  on  the  computer.  Evaluation  in  several 
typical  cases  showed  that  the  variance  of  the  estimate  using  the  sin  x/x  interpolator 
was  2  to  3  times  that  of  the  quadratic  Interpolator,  in  spite  of  the  fact  that  the 
quadratic  interpolator  is  much  simpler.  Both  interpolators  were  used  to  process 
sea  tape  data  described  In  Section  4.0,  and  the  results  of  this  processing  supports 
the  analysis. 

3.5  Cancellation  of  Interferences  Prior  to  Tracking.  There  are  basically  two 
approaches  to  dealing  with  interference  in  the  adaptive  tracking  system.  One  is  to 
depend  upon  the  spatial  response  of  the  split  arrays  and  upon  the  adaptive  proper¬ 
ties  of  the  tracker  to  provide  relative  immunity  to  the  interference.  The  other  ap¬ 
proach  is  to  attempt  to  cancel  the  interference  from  the  half  beam  inputs  to  the 


adaptive  tracker.  It  was  shown  during  phase  2  of  this  study  C  23  that  the  primary 
effect  of  an  interference  is  to  introduce  a  bias  in  the  bearing  estimate,  and  that 
this  bias,  can  be  quite  severe  when  the  signal-to-interference  ratio  is  low  at  the 
split  array  outputs.  Under  these  circumstances,  it  may  be  desirable  to  cancel  the 
interference  prior  to  tracking. 

The  adaptive  noise  canceller  C  53  is  often  used  to  eliminate  an  interference 
arriving  on  the  sidelobes  from  a  beam  output.  A  reference  sensor,  for  example 
an  omnidirectional  hydrophone  that  is  spatially  separated  from  the  array,  is  used 
as  a  reference  in  the  LMS  canceller  configuration  of  Figure  11.  When  the  interfer¬ 
ence  power  dominates,  the  adaptive  filter  will  spatially  reject  the  interference 
from  the  beam  output.  Depending  upon  the  spacing  between  the  reference  hydro¬ 
phone  and  array  and  upon  the  separation  of  the  target  and  interference,  some  re¬ 
jection  of  the  target  signal  will  occur,  but  this  will  usually  be  acceptable  in  light 
of  the  interference  rejection.  Figure  12  shows  this  approach  applied  to  the  adaptive 
tracker.  A  reference  is  supplied  for  each  half  array  and  a  canceller  implemented 
in  each  half  array  output.  A  convenient  choice  for  a  reference  may  be  a  single 
hydrophone  located  in  the  opposite  half  array. 

Appendix  F  considers  the  possible  bias  of  the  broadband  and  narrowband  bear¬ 
ing  estimates  due  to  the  transfer  functions  of  the  adaptive  cancellers  prior  to  the 
tracker  in  Figure  12.  This  is  done  by  determining  the  steady  state  mean  weight  of 
the  adaptive  trader  via  a  Wiener  filter  approach.  A  bias  in  the  broadband  estimate 
is  indicated  by  a  shift  in  the  position  of  the  peak  of  the  mean  weights,  while  a  bias 
in  the  narrowband  estimate  appears  as  a  bias  in  the  phase  of  the  frequency  domain 
weight  containing  the  signal.  It  was  assumed  in  all  cases  that  the  Interference  power 
dominated  the  signal  and  background  noise,  since  this  is  the  environment  in  which 
the  canceller  is  applied. 


In  the  broadband  case,  Appendix  F  shows  that  when  the  signal  to  background 
noise  ratio  is  large  (both  dominated  by  interference),  the  peak  in  the  weights  due 
to  the  interference  is  effectively  cancelled  with  no  bias  in  the  signal  peak,  hence, 
no  bias  in  the  broadband  bearing  estimate.  The  amplitude  of  the  signal  peak  after 
cancelling  is  essentially  that  of  the  Interference  peak  prior  to  cancelling,  and  is 
therefore  easily  tracked.  When  the  signal  to  background  noise  ratio  is  small,  the 
interference  is  still  effectively  rejected,  but  a  small  bias  is  introduced  into  the  sig¬ 
nal  peak.  This  results  in  a  delay  estimate  bias  of 

®T  “  '  ^  (*  *  I)  |Sin  9S  ~  SiD  ®l|  (35) 

where 

Bj  =  split  beam  response  in  direction  of  the  interference  (sidelobe  response) 
Bg  =  split  beam  response  in  direction  of  the  signal 
D  =  distance  between  split  array  phase  centers 
d  =  distance  from  split  array  to  its  reference 

*  signal  bearing 

*  Interference  bearing. 

C  *  speed  of  sound 

Since  the  sidelobe  response  is  small  in  comparison  to  the  malnlobe  response, 

Bg  >>  Bj,  B^  Is  small.  It  can  be  made  arbitrarily  small  by  placing  the  references 
near  the  array  phase  centers  so  that  d/D  **  1.  In  addition  to  the  bias  introduced  in 
the  low  SNR  case  the  entire  weight  vector  is  attenuated  proportionally  to  the  SNB 
at  the  beamformer  output,  so  longer  time  constants  will  be  required  if  degradation 
in  estimation  variance  is  not  tolerable. 

In  the  narrowband  case,  tracking  is  only  possible  at  relatively  high  SNR,  as 


shown  in  ( 2].  Under  this  condition,  Appendix  F  shows  that  the  bias  In  the  phase 
of  the  frequency  domain  weight  containing  the  signal  is  approximately 
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with  Pg  and  Pj  the  signal  and  interference  powers,  respectively.  Since  both  Pg/Pj 
and  Bj/Bg  are  small,  the  residual  bias  in  the  narrowband  estimate  is  small. 

3.6  Simulations  of  an  Adaptive  and  a  Conventional  Split  Beam  Tracker  With  a 
Moving  Target.  During  the  first  phase  of  the  adaptive  tracker  study  ^  the 
adaptive  tracker  performance  was  compared  to  a  conventional  adaptive  tracker 
described  in  C6j  with  static  targets  Figure  13  shows  a  block  diagram  of  the 
tracker.  The  input  signals  from  the  left  and  right  half  beams  are  hard  clipped  and 
sampled,  with  the  clipped  signal  for  one  side  subjected  to  a  variable  delay,  r.  This 
delayed  signal  and  the  other  clipped  input  are  processed  by  a  two  point  correlator, 
as  shown,  which  computes  the  correlation,  0(t),  of  the  clipped  signal  at  r  -  T. 

The  value  of  0(t)  is  either  +1  or  -1  for  the  clipped  inputs.  The  contents  of  the  de¬ 
lay  register,  which  determines  the  value  of  the  adjustable  delay,  is  increased  by 
At  when  +1  occurs  on  input  2  and  -1  on  input  1.  It  is  decreased  by  At  when  +1 
occurs  on  input  1  and  -1  on  input  2,  with  no  change  in  count  when  the  two  agree. 

The  tracker  is  In  steady  state  when  4(t  -  T)  and  0(t  +  T)  are  equal  (both  +1  or 
both  -1).  Assuming  that  the  correlation  function  of  the  input  is  symmetrical,  this 
means  that  the  value  of  t  is  equal  to  the  delay  between  the  two  split  beam  inputs. 

The  analytical  comparison  of  the  two  trackers  using  the  results  of  C  l]  and 
C 63  indicated  that  the  variance  of  the  adaptive  tracker  bearing  estimate  was  nom¬ 
inally  a  factor  of  3  improved  over  the  conventional  tracker.  This  was  supported 
by  computer  simulations  of  the  two  trackers.  This  appendix  repeats  these  simula¬ 
tions  with  a  dynamic  target.  The  signal  has  linearly  increasing  time  delay  between 
array  halves,  which  is  the  same  as  linearly  varying  bearing  for  a  target  at 
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broadside.  Both  the  target  and  the  background  noise  are  uncorrelated  gaussian 
sequences,  and  the  background  noise  between  the  array  halves  is  also  uncorre¬ 
lated.  The  algorithm  sample  rate  is  24,000  Hz  and  no  interpolation  is  used  be¬ 
tween  taps  of  the  adaptive  filter. 

The  statistics  of  the  bearing  estimate  were  computed  from  an  ensemble  of 
10  runs,  with  time  averaging  performed  over  4. 67  seconds  in  each  run  of  the 
ensemble.  The  procedure  was  as  follows.  For  a  given  signal-to-noise  ratio,  a 
value  of  the  feedback  coefficient,  p,  was  chosen  so  that  the  adaptive  tracker  could 
acquire  and  track  the  signal  throughout  the  duration  of  the  simulation.  Since  the 
goal  was  to  compare  the  trackers  during  tracking,  any  runs  with  spurious  loss  of 
track  were  not  used.  The  4. 7  second  time  averages  of  estimate  bias  and  variance 
were  computed  for  each  such  run,  and  an  ensemble  of  ten  of  these  runs  averaged 
to  yield  a  value  of  bias  and  variance  for  a  given  signal-to-noise  ratio  and  bearing 
rate. 

In  order  to  compare  the  two-point  correlator  tracker  of  C6D  to  the  adaptive 
tracker,  its  feedback  parameter,  A,  was  varied  experimentally  until  the  bias  of 
the  two-point  correlator  estimate  was  approximately  the  same  as  that  of  the  adap¬ 
tive  tracker  for  the  given  signal-to-noise  ratio  and  bearing  rate.  As  in  the  case  of 
the  adaptive  tracker,  all  statistics  for  the  two  point  correlator  were  based  upon 
an  ensemble  average  of  simulation  runs,  with  time  averaging  over  4.67  seconds  in 
each  run.  Due  to  greater  run-to-run  fluctuations  in  the  estimate  from  the  two  point 
correlator,  twenty  runs  were  included  in  an  ensemble.  With  the  bias  of  the  two 
trackers  set  approximately  equal,  it  is  possible  to  compare  the  performance  based 
upon  bearing  estimate  variance  at  a  given  signal-to-noise  ratio,  bearing  rate,  and 
bias  (or  lag). 
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Table  1  shows  the  results  of  these  simulations  for  a  signal-to-aoise  ratio  of 
0  dB  and  three  rates  of  change  in  delay  between  array  halves.  Table  2  converts 
these  delay  rates  to  bearing  rates  for  two  arrays,  one  with  a  7. 5  foot  baseline  and 


the  other  with  a  75  foot  baseline.  The  feedback  coefficient  of  the  adaptive  filter,  , 
—18 

has  been  set  to  2  to  assure  tracking.  The  standard  deviation  of  the  adaptive 

tracker  estimate  ranges  from  0.751  to  0.975  times  that  of  the  two-point  correlator. 

Table  3  shows  the  same  statistics  for  a  -10  dB  slgnal-to-noise  ratio,  with 
-18 

H  ■  2  .At  this  lower  signal  to  noise  ratio,  the  advantage  of  the  adaptive  tracker 
is  greater,  with  its  standard  deviation  ranging  from  0. 535  to  0. 73  times  that  of 
the  two  point  correlator. 

It  is  interesting  to  note  that  the  advantage  demonstrated  by  the  adaptive 
tracker  in  these  simulations  is  quite  close  to  that  predicted  analytically  in 


TABLE  1.  COMPARISON  OF  ADAPTIVE  TRACKER  AND  TWO  POINT 
CORRELATOR,  SIMULATIONS  WITH  DYNAMIC  TARGET:  SNR  =  0  dB,  H  -  2~16 


Rate 

Bias 

Standard  Deviation 

Ratio 

M  sec/sec 

Adaptive 

Tracker 

Two  Point 
Correlator 

Adaptive 
Tracker,  oadt 

Two  Point 
Correlator, 

<radt//<rtpc 

65.45 

6.508  x  10-5 

6.6  x  10"5 

2.369  x  10"5 

2.429  x  10-5 

0. 975 

130.9 

7.846  x  10-5 

7.900 x  10"5 

2.394  x  10"5 

2.804  x  10'5 

0.854 

261.8 

9.633  x  10-5 

9.241 x  10"5 

2.393  x  10"5 

3. 184  x  10"5 

0. 751 

TABLE  2.  CONVERSION  OF  DELAY  RATES  TO  BEARING  RATES  FOR  TWO 

HYPOTHETICAL  ARRAYS 


Delay  Rate 

Bearing  Rate 

7.5  ft  Baseline 

Bearing  Rate 

75  ft  Baseline 

65.45  s/s 

2. 5° /sec 

0.25°/sec 

130.9  s/s 

5° /sec 

0.5° /sec 

261. 8  s/s 

10° /sec 

l°/sec 
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TABLE  3.  COMPARISON  OF  ADAPTIVE  TRACKER  AND  TWO  POINT 
CORRELATOR,  SIMULATIONS  WITH  DYNAMIC  TARGET:  SNR  =  10  dB,  n  *  2“  18 


Rate 

Bias 

Standard  Deviation 

Radio 

H sec/sec 

Adaptive  Two  Point 

Tracker  Correlator 

Adaptive 
Tracker,  <x&dt 

Two  Point 
Correlator, 

^adt^tpc 

65.45 

8.365  x  10"5  8.023  x  10”5 

4. 590  x  10“5 

8.576  x  10"5 

0. 535 

130.9 

1.239 x  10-4  1.265  x  10"4 

8.03  x  10"5 

1.422  x  10*4 

0*  o65 

261.8 

1. 989  x  10~4  1.954  x  10“4 

1  ■  ,  - 

8.722  x  10-5 

1. 194  x  10"4 

0.730 

Appendix  G  of  C2j.  For  a  somewhat  different  split  beam  tracker,  the  adaptive 
tracker  standard  deviation  was  predicted  to  be  0. 831  times  that  of  the  split  beam 
tracker  under  the  condition  that 
„  2 

be  4  «  |i<r  n  <*:  1 
where 

b  =  input  bandwidth 
c  3  rate  of  change  in  delay 
<5  =  algorithm  sample  rate 
2 

<7q  *  input  power. 

2  2 

For  these  simulations,  par n  1  and  (bcfi/^O’jj)  varies  from  0.214  to  3.436.  Since 

the  referenced  analysis  used  a  bearing  deviation  indicator  structure  for  the  con¬ 
ventional  tracker,  and  an  exponential  correlation  function  for  the  inputs,  these 
simulations  should  be  used  only  for  qualitative  validation  of  the  parametric  trends 
in  that  analysis. 


SECTION  4.  0 

ADAPTIVE  TRACKER  PERFORMANCE  WITH  BROADBAND 
SEA  TAPE  DATA 


4.0  ADAPTIVE  TRACKER  PERFORMANCE  WITH  BROADBAND  SEA  TAPE  DATA 

As  part  of  the  third  phase  of  the  adaptive  tracker  study ,  several  tapes 
recorded  at  sea  using  a  towed  line  array  were  processed  using  the  broadband  adap¬ 
tive  tracker  algorithm.  These  sea  tests  utilized  broadband  random  noise  sources 
suspended  from  one  or  two  drifting  surface  vessels,  with  the  array  moving 
relative  to  the  sources.  Two  such  encounters  were  processed,  a  "HOTEL" 
tape  in  which  only  a  single  source  was  present ,  and  an  "INDIA"  tape  which 
included  two  sources  simultaneously.  The  sources  are  extremely  strong, 
dominating  background  and  tow  ship  noise. 

The  recording  consists  of  17  hydrophone  cluster  outputs  from  the  array, 
with  the  clusters  spaced  every  7.25  feet.  The  hydrophone  outputs  are  band- 
limited  to  500  Hz  and  sample  at  1  KHz.  Figure  14  shows  an  estimate  of  the 
power  spectrum  of  the  output  of  a  single  hydrophone  during  the  "HOTEL" 
exercise,  averaged  over  32.7  seconds.  Figure  15  shows  the  corresponding 
estimate  of  the  autocorrelation  function  (ACF)  of  the  hydrophone  output. 

The  main  peak  is  only  several  taps  wide  (due  to  sampling  at  the  Nyquist  rate 
with  a  very  sharp  bandlimiting  filter)  and  two  large  negative  peaks  appear 
adjacent  to  the  main  peak.  The  shape  of  the  ACF  is  important  in  that  the 
weight  vector  of  the  tracker  will  have  essentially  the  same  shape,  and  the 
bearing  estimate  is  obtained  by  locating  the  peak  of  the  weights. 

Because  of  the  high  signal  to  noise  ratio  during  the  exercises,  and  the 
relatively  low  directivity  of  the  half  arrays  over  much  of  the  band ,  single 
hydrophones  were  used  as  the  inputs  to  the  adaptive  tracker  rather  than 
split  arrays.  The  hydrophones  selected  provide  a  baseline  of  72.5  feet. 

Two  types  of  runs  were  made  with  the  "HOTEL"  data,  processing  the  raw 
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Figure  15.  Estimated  Autocorrelation  Function  of 
Hydrophone  #4  Output,  Exercise  HOTEL- 1 


hydrophones  outputs  first ,  then  tracking  with  additional  computer  generated 
noise  added  to  evaluate  the  tracker  at  lower  signal- to- noise  ratios.  The 
"INDIA”  encounter  was  similarly  processed. 

4.1  Processing  of  Unmodified  "HOTEL"  Data.  The  geometry  of  the 
"HOTEL-1"  exercise  is  shown  in  Figure  16.  The  array  is  towed  past  the 
noise  source  at  2.8  knots  at  a  depth  of  350  feet.  The  source  is  suspended 
at  100  feet  with  a  nominal  spectrum  level  of  160  dB  //y  Pa  at  one  yard.  The 
closest  point  of  approach  is  approximately  2214  yards,  as  shown  in  Figure  16. 
Assuming  cylindrical  spreading  and  a  sea  state  of  between  2  and  3,  this  would 
suggest  an  SNR  of  20  to  30  dB  at  the  hydrophones  in  the  120-480  Hz  band. 

The  bearing  estimate  vs.  time  for  three  segments  of  the  HOTEL- 1  tape 
are  shown  in  Figures  16  through  19.  The  feedback  coefficient,  y ,  has  been 
set  at  2~10,  since  the  signal-to-noise  ratio  is  high.  The  variance  of  the 
estimate  is  clearly  weil  under  0.1  degrees  throughout  the  exercise,  which  is 
not  unreasonable  given  the  high  signal-to-noise  ratio. 

Figure  21  shows  the  weight  vector  vs.  time  for  a  64  tap  adaptive  tracker 
with  y  =  2-10  during  the  period  6:39  to  6:46  of  the  HOTEL-1  exercise. 

The  track  of  the  noise  source  is  clearly  visible.  Table  4  compares  the  esti¬ 
mated  bearing  of  the  source  from  the  adaptive  tracker  with  the  bearing 
reconstructed  graphically  from  Figure  16.  The  "HOTEL- 1"  tape  has  also 
been  used  to  evaluate  several  beamforming  schemes,  [7] .  Using  the  result 
of  these  tests,  given  as  a  BTR  plot  in  Figure  20,  the  apparent  bearing  of 
the  target  can  be  seen  to  be  between  about  33°  and  20°  over  the  same  period 
of  time,  consistent  with  the  adaptive  tracker  re  stilts.  Part  of  this  difference 
between  apparent  and  reconstructed  bearings  will  be  due  to  the  difference 
in  depth  between  array  and  source .  Reference  [  7]  has  also  observed  a  yaw 
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’  Figure  18.  Bearing  Estimate  from  Adaptive  Tracker  for  HOTEL  Run 
Start  of  run  *  06:43:18.41 


Figure  20.  Bearing  Time  Plot  for  HOTEL- 1  Exercise  with  Beamfbrming 
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angle  of  5°  -  10°  during  the  test ,  and  indicated  the  possibility  of  array  pitch 
and  deviation  of  array  depth  from  the  desired  value.  Further,  local  varia¬ 
tions  in  array  shape,  such  as  array  curvature,  could  cause  bearing  errors, 
particularly  when  only  two  hydrophones  are  used  in  the  bearing  estimate. 

By  using  a  different  pair  of  hydrophones  with  the  same  spacing  as  the  input 
to  the  adaptive  tracker  with  the  identical  input  data,  it  was  shown  in  the  . . 
second  quarter  that  array  curvature  does  not  seem  to  be  a  component  in  the 
bearing  error,  since  the  bearing  estimates  were  virtually  unchanged.  This 
conclusion  is  reinforced  by  the  agreement  between  the  adaptive  tracker 
results  and  Figure  17,  where  the  beamforming  [7]  used  all  hydrophones,  and 
would  be  relatively  immune  to  local  variations  in  array  shape. 

If  a  fixed  bias  of  about  11°  is  subtracted  from  the  estimates,  as  shown 
in  Table  4,  excellent  agreement  is  obtained  with  the  reconstructed  values. 
This  strongly  suggests  that  t  bias  of  approximately  11°  is  the  combined 
effect  of  the  above  factors. 

4.2  Processing  of  HOTEL  Data  at  Reduced  Signal-to-Noise  Ratio.  In 
order  to  use  the  "HOTEL- 1"  data  at  lower  signal-to-noise  ratios ,  computer 
generated  Gaussian  noise  was  added  to  the  two  hydrophone  outputs.  The 
noise  is  uncorrelated  in  time  and  between  the  two  inputs.  Since  the  exact 
signal-to-noise  ratio  on  the  tape  is  unknown ,  the  power  ratio 

P„  +  p„ 

R_  =  %  ■■  n  (66) 

P  PA 

was  varied,  where  P.  is  the  signal  power  on  the  tape,  P_  the  noise  power 
on  the  tape,  and  P^  the  additional  noise  power  added.  The  actual  signal- 


to-noise  ratio  is  then 


with  (P  /P  )  the  unknown  signal-to- noise  ratio  on  the  tape, 
s  n 

This  modified  acoustic  data  was  applied  to  the  adaptive  tracker  structure, 
a  train  a  64  tap  time  domain  filter,  and  v  varied  from  2-10  to  2-16.  In  addi¬ 
tion,  the  two  interpolators  analyzed  in  Appendix  E  were  used  so  that  their 
performance  could  be  compared  with  real  data.  The  results  were  averaged 
over  26.6  seconds  to  produce  a  sample  mean  and  variance  for  the  bearing 

estimate.  Table  5  shows  the  mean  bearing  for  7  values  of  R  and  for  values 

P 

of  u  for  both  the  quadratic  and  sin  x/x  interpolator,  while  Table  6  shows 
the  variance  for  the  same  runs.  The  average  is  over  the  period  of  6:44:38 
to  6:45:04  of  the  exercise  during  which  the  reconstructed  bearing  will  be 
nominally  14°.  These  results  show  generally  good  agreement  with  the 
theoretical  predictions  for  a  static  broadband  target  developed  in  Section  II . 
Some  deviation  from  these  predictions  is  to  be  expected  since  the  target  has 
a  slow  bearing  rate  and  the  spectrum  is  not  exactly  white.  Also,  these  runs 
indicate  that  the  sin  x/x  interpolator  has  approximately  1.5  times  the 
variance  of  the  quadratic  interpolator,  which  is  in  the  theoretically 
predicted  range. 


e  5.  Standard  Deviation  of  Bearing  Estimate  for  HOTEL- 1  Tape 


4. 3  Processing  of  the  INDIA  Data.  In  the  INDIA  exercise,  the  array  is  towed 
between  two  suspended  broadband  sources,  in  the  geometry  shown  in  Figure  22. 
Projector  1  in  the  figure  is  radiating  broadband  noise  at  a  nominal  source  level 
of  160  dB  re  1  pPa  at  1  yard  over  the  17  to  500  Hz  band,  addle  projector  number 
2  ie  operating  at  a  source  level  of  120  dB.  The  stronger  source  range  in  nomin¬ 
ally  J  Kyds,  while  the  weaker  source  is  a  nominally  3  Kyds,  so  the  projector  1 
will  appear  nominally  36  dB  stronger  at  the  array.  The  geometry  is  such  that 
projector  1  appears  initially  at  60°  forward  of  broadside  and  is  at  48°  forward  of 
broadside  at  die  end  of  the  exercise.  Projector  2  starts  at  32°  forward  of  broad¬ 
side  and  finishes  14°  forward.  The  portion  of  Ihe  exercise  processed  by  the 
adaptive  tracking  extends  from  0035  to  0052. 

The  outputs  of  hydrophones  4  and  14  were  applied  to  a  broadband  adaptive  filter 
structure  with  64  taps  and  with  y  set  at  2~10.  The  adaptive  filter  weights  vs  time, 
plotted  in  Figure  23,  show  three  primary  peaks.  The  largest  starts  at  about  60° 

(or  tap  45),  with  a  smaller  peak  slightly  to  the  right  moving  parallel  to  the  large 
peak.  Reference  [7]  notes  that  the  peak  associated  with  the  strong  source  is  often 
a  doublet,  and  that  a  peak  of  unknown  origin  appears  approximately  15°  to  the 
right  of  the  TPafo  peak,  moving  parallel  to  it.  Ihese  are  possibly  surface  reflec¬ 
tions  of  projector  1.  The  third  peak,  considerably  weaker  than  the  main  peak 
starts  at  about  28°,  or  tap  39.  This  appears  to  be  projector  2,  based  upon  re¬ 
construction  from  Figure  22  and  the  information  provided  [7]. 
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Figure  22.  Array /Source  Geometry  (or  India  Exercise 
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Figure  23.  Adaptive  Tracker  Weights  vs  Time  for  INDIA  Data 


The  adaptive  tracker  was  assigned  to  the  peaks  at  tap  45  and  tap  39 ,  using 
the  windowing  scheme  described  in  Section  2. 3.  The  window  was  selected  to  be 
3  taps  wide,  centered  on  the  tap  nearest  the  previous  delay  estimate.  However, 
the  window  center  was  not  allowed  to  move  more  than  one  tap  per  iteration  to  pre¬ 
vent  loss  of  track  due  to  fluctuations.  Figures  24  and  25  show  the  trades  developed 
for  the  two  sources  projectors  1  and  2,  respectively.  Table  5  compares  the 
estimates  with  the  bearing  reconstructed  from  the  tracker  of  22.  The  estimates 
for  projector  1  tend  to  be  high  and  relatively  constant,  although  the  reconstructed 
bearing  is  decreasing.  The  estimates  are,  however,  consistent  with  the  results 
of  [7],  in  which  the  peak  in  the  beamformer  response  due  to  projector  1  was  con¬ 
stant  at  approximately  59°.  This  indicates  that  some  external  phenomenon,  such 
as  array  geometry,  multipath,  or  interference  between  the  two  targets  accounts 
for  the  shift  in  the  estimates.  Reference  [7]  notes  that  during  this  period  a  prom- 
inant  surface  reflection  from  the  tow  ship  appears  at  60°,  so  it  is  assumed  that 
this  is  the  cause  of  die  bias.  The  track  for  projector  2  shows  excellent  agreement 
with  the  reconstructed  bearings. 

As  in  section  4.2,  the  INDIA  tapes  were  processed  with  additional  computer 
generated  random  noise  added  to  die  hydrophone  output.  Since  the  signal-to-noise 
ratio  on  the  tape  was  unknown,  the  power  ratio,  Rp  given  by  (66)  was  varied  from 
-10  dB  to  +10  dB  in  10  dB  steps.  An  Rp  of  X  dB  means  that  the  total  signal  plus 
noise  power  on  the  tape  was  X  dB  above  the  added  noise.  Figures  26  and  27  show 
die  tracks  with  Rp  -  10  dB,  Figures  28  and  29  with  Rp  *  0  dB,  and  Figures  30  and 
31  with  Rp  a  -io  dB.  The  tracks  are  maintained  at  Rp  =  10  dB  and  0  dB  for  both 
contacts.  However,  when  Rp  ■  -10  dB,  the  tracks  are  eventually  lost,  even 
though  the  contacts  appear  to  be  maintained  initially. 
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Figure  24.  Track  (or  Projector  1,  India  Data  Start  of  Run  ■  0035,  No  Added  Noise 


Figure  25.  Track  for  Projector  2,  India  Data  Start  of  Run  *  0035,  No  Added  Noise 
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Figure  26.  Track  for  Projector  1,  India  Data  Start  of  Run  ■  0035,  Computer  Generator  Noise 
Added,  Rp  -  lOdB 


07010-33 
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Figure  30.  Track  for  Projector  1,  India  Data  Start  of  Run  »  0035,  Computer  Generated 
Noiae  Added,  Rp  *  -10  dB 
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SECTION  5. 0 


CONCLUSIONS 


5.0  CONCLUSIONS  -  PHASES  1,  2  AND  3 

A  new  split  beam  tracker  structure  has  been  presented,  analyzed  and 
simulated  under  a  variety  of  single  target,  multiple  target,  broadband, 
narrowband,  static  and  dynamic  situations.  Phase  1  demonstrated  the 
theoretical  feasibility  of  the  device  under  idealistic  input  conditions ,  showing 
close  agreement  with  lower  bounds  on  the  variance  of  the  bearing  estimate. 

In  Phase  2,  the  performance  in  more  realistic  environments  was  analyzed 
and  simulated.  In  particular,  target  dynamics  and  changing  input  spectra 
that  differ  from  those  used  in  the  design  of  a  conventional  tracker  were 
studied.  The  tracker  responds  much  faster  than  the  convergence  time  of  the 
adaptive  filter  itself  since  the  largest  of  the  weights  can  be  selected  far  before 
the  weights  have  achieved  steady-state.  In  the  comparisons  using  changing 
input  spectra,  the  adaptive  filter  tracker  demonstrated  the  potential  for  drama¬ 
tic  improvements  over  a  conventional  tracker  which  is  operated  in  an  environ¬ 
ment  other  than  for  which  it  was  designed.  This  phase  concluded  that  the 
adaptive  filter  tracker  has  all  of  the  performance  qualities  one  would  associate 
with  a  tactical  broadband  or  narrowband  tracker,  thus  logically  leading  to  the 
third  phase. 

Phase  3  addressed  the  systems  considerations  in  the  application  of  the 
adaptive  tracker  to  sonar  arrays .  This  included  the  simultaneous  tracking  of 
multiple  targets,  the  effects  of  strong  plane  wave  interference,  the  combining 
of  both  broadband  and  narrowband  target  data  ( rather  than  separating  the 
data  into  a  broadband  tracker  and  a  narrowband  tracker) ,  the  effects  of 
arrays  on  the  inputs,  and  general  operational  system  considerations.  Again, 
the  analytical  results  were  verified  with  computer  simulations.  In  addition, 
recorded  sea-data  from  the  elements  of  a  towed  array  were  employed  as  well. 
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In  summary,  there  are  some  genuine  advantages  in  performance  and 
perhaps  even  in  implementation  offered  by  the  adaptive  niter  tracker  over 
conventional  trackers .  The  performance  benefits  come  when  signal  or  noise 
statistics  are  not  known  to  the  designer  of  the  front-end  filters  on  a  conven¬ 
tional  system.  Implementation  benefits  can  occur  when  sampling  rates  at  the 
beamformer  outputs  are  much  slower  than  the  time  resolution  required  to 
achieve  a  specified  angular  estimation  accuracy.  >  The  adaptive  tracker  inter¬ 
polates  between  tap  values  rather  than  having  to  sample  at  a  higher  rate.  ^ 
In  terms  of  dynamic  behavior,  the  adaptive  filter  tracker  is  understood, 
analysis  predicting  simulated  tracking  rates  and  the  resulting  possible  lag. 
Windowing  on  the  filter  weights  allows  handling  of  multiple  targets,  and 
putting  second-order  dynamics  on  the  track  window  functions  enables  a 
single  adaptive  filter  tracker  to  follow  more  than  one  target,  even  if  their 
tracks  cross.  Although  the  signal-to-noise  ratio  performance  is  comparable 
to  that  of  nonadaptive  trackers  under  ideal  conditions,  in  practical  unknown 
environments  the  adaptive  filter  tracker  can  realize  the  predicted  perform¬ 
ance  via  its  capabilities  to  estimate  statistics  and  adjust  accordingly . 

Thus,  in  the  three  phases  of  this  study,  the  LMS  adaptive  filter  con¬ 


figured  as  a  bearing  tracker  has  been  carried  from  a  concept ,  to  a  well 
understood  design,  to  a  stage  ready  for  system  application,  test  and 
evaluation. 
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APPENDIX  A 


EFFECTS  OF  MULTIPATH 


In  the  work  during  the  first  two  phases  of  the  adaptive  tracking  systems  study 
[1,2],  the  signal  has  been  assumed  to  arrive  at  the  sonar  array  via  a  single  plane 
wavefront.  In  most  sonar  environments,  however,  the  signal  arrives  via  multiple 
paths,  for  example,  direct,  surface  reflection,  and  bottom  bounce.  This  appen¬ 
dix  considers  the  effects  of  the  multipath  propagation  on  the  performance  of  the 
split  beam  adaptive  tracker  and,  in  particular,  assesses  the  sensitivity  of  the 
adaptive  tracker  to  multipath  relative  to  conventional  trackers. 

A  model  for  the  outputs  of  the  split  array  beamformers  with  N  multipath 
arrivals  of  a  single  target  is 
N 

X(t)  =  5j  a  S(t-Y )  +  n  (t)  (A-l) 

n=l  n  n  x 

and 


N 

d(t)  =  £  bns(t"Xn'Tn)  +  nd(t)  (A'2) 


Here,  Sq  andbn  are  the  attenuations  of  the  signal  along  the  nth  path  to  the  left 
half  and  right  half  arrays ,  respectively,  VQ  and  \n  are  the  propagation  delays 
along  those  paths,  and  rQ  is  the  inter-array  delay  associated  with  each  path.  The 
value  of  rn  depends  upon  the  angle,  )3q,  the  wavefront  makes  with  the  axis  of  the 
phase  centers. 


(A-3) 


where 


sin0Q  =  cos$nstn0  (A-4) 

with  6  the  target  bearing  for  all  paths  and  0Q  the  arrival  angle  of  the  nth  path  in 
the  vertical  plane  at  angle  6 . 

This  problem  can  be  looked  at  in  terms  of  linear  system  theory  as 

x(t)  =  hx(t)^(t)  +  nx(t)  (A -5) 


and 

d(t)  =  hd(ty%(t)  +  nd(t)  (A-6) 

where  x%  is  the  convolution  of  x  and  y  and  where  hx  and  hd  are  impulse  responses 
given  by 

Vrt-r  ssc-y  (A-T) 

N 

=  £  b  8<T-X  -T  )  (A-8) 

n=l 

Let  H  (w)  and  H  ,(w)  be  the  transfer  functions  corresponding  to  h  (r)  and  h  ,(t). 

X  0  X  Cl 

Then  the  power-spectral  density  of  x(t)  is 

»„<“»  *  sra(">  (A-9) 

and  the  cross-spectrum  between  d  and  xls 

Sdx*(w)  =  Hd(w)Hi(w)  Sss(w)  {A"10) 

with 

S  (w)  =  power  spectrum  of  s(t) 

S^fw)  »  power  spectrum  of  nx(t) 


The  adaptive  tracker  weight  will  converge  to  the  discrete  Wiener  filter  for 


these  inputs,  which  approximates  the  continuous  Wiener  filter  if  the  filter  is 
long  enough  and  the  sample  rate  Is  adequate.  The  continuous  Wiener  filter  for 
this  example  Is  given  by  the  ratio  of  (A -10)  to  (A -9) 


Wf(w) 


Hd(fa/)H*(t>?)Saa(h)) 
|hv(w)|2S,s(o>)  +  S^w) 


(A-ll) 


X'-'1  ss 
From  (A-7)  and  (A-8) 

H  («)  =  Z  a  e'j“Yn 

x  n=l  “ 


(A-12) 


and 

Hj(w)  -  £  h.-*"  W  «A-13» 

n=l  n 

In  most  sonar  encounters,  the  dimensions  of  the  array  are  small  in  com¬ 
parison  to  the  range  to  the  target,  so  that  the  target  is  in  the  far  field  (otherwise 
the  plane  wave  assumption  is  invalid).  It  is  therefore  reasonable  to  assume  that 
the  array  is  a  point  receiver  with  respect  to  the  target,  so  that  the  propagation 
delays  along  any  path  to  the  half  arrays  are  the  same,  i.e. ,  Yn=A.Q.  Consequently, 

a_=b_  and 
n  n 

N 

H  (w)  =  Z  a  e'jWVn  (A-14) 

35  n=l  n 


and 


Hd(«)  = 


E  »  e-l"<vlc‘,,V‘n8' 

ftl  “ 


(A- 15) 


One  comparison  of  interest  is  between  the  adaptive  tracker  and  a  conventional 
split  beam  tracker  which  bases  its  delay  estimate  on  the  location  of  the  peak  of 
the  correlation  function  of  the  split  beam  outputs  (or,  equivalently,  on  the  zero 
crossing  of  its  derivative >[4] .  If  Hw(u)  is  the  transfer  function  of  the  fixed 
spectral  weighting  filter  of  the  conventional  tracker,  the  cross 'spectral  density 
between  the  split  array  outputs  is 


C(«)  =  |Hw(o>)|2  Hd(w)H*(w)Sga(a>) 


(A- 16) 


From[4],  the  optimal  choice  of  |H  (w)|  is  the  generalized  Eckart  filter, 

|hwM2  - - ^ - 


SQ(o;)+MSgs(w)Sn(w) 


(A- 17) 


Here  it  has  been  assumed  that  the  details  of  the  multipath  structure  are  unknown 
and  the  Eckart  filter  design  must  be  based  upon  the  signal  and  noise  spectra 
only.  At  low  signal-to-noise  ratio  (SNR),  (A-ll)  reduces  to 


Wf(aj>~S— M  Wss<w» 


(A-18) 


nx 


while  (A-16)  becomes 


Sss(w) 

CM  =‘-2| —  H.(«>H  *(u)S  («) 
g2  d  x  ss 


(A-19) 


Now,  assume  that  the  signal  and  noise  both  have  the  same  spectral  shape,  i.e. , 

Sg(a>)  =  P§SM  and  SQ(w)  =  PQS(«)  (A-20) 

so 


Wf(w)*  p  HdMHxM 

n 


(A-21) 
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and 


P2 

C(cj)  =  —  Hd(oj)H  *M 
Pa 


(A-22) 


Comparing  (A-21)  and  (A-22),  it  can  be  seen  that  at  low  SNR,  the  adaptive  filter 
weight  vector  and  the  cross-correlation  function  have  the  same  form  except  for  a 
scale  factor.  Therefore,  the  two  trackers  will  exhibit  the  same  sensitivities  to 
multipath,  at  least  in  terms  of  mean  behavior. 

Similarly,  a  conventional  narrowband  estimate  would  be  based  upon  the  phase 
of  the  cross-spectrum  at  the  signal  frequency,  uQ.  where  high  SNR  is  of  interest 
for  most  broadband  targets.  From  (A-ll)  and  (A-16),  it  can  be  seen  that  the 
weight  vector  and  cross -spectrum  will  have  the  same  phase,  so  that  the  sensi¬ 
tivity  of  the  two  trackers  to  multipath  will  be  the  same. 

In  order  to  examine  the  sensitivity  of  the  adaptive  tracker  to  multipath  in 
the  broadband  case  (low  SNR),  let 


Sas(W> 

and  let  g(r)  be  the  Inverse  Fourier  transform  of  G(u).  Also,  let 

N  N 

Vw)  -  Hd(a»H*(w)  -  £  2  aneme 

n=l  m=l 


(A-23) 


-  Tn)  +  Ta 


(A-24) 


and  its  inverse  transform  be 


N  N 

VT)  "  2  ^  Vmi(T "  fym  "  V  "  V  <A-25) 

n=l  m=l 

Using  (A-18),  the  adaptive  filter  weight  vector  Is  approximately  given  by  the 
convolution  of  g(r)  and  k^fr),  which  is  a  weighted  sum  of  shifted  versions  of  g(r), 
i«  e#  i 

A- 


(A-26) 


Wf(T)  =  Z  Z  anV(T  -  (\  "  V  -  Tn> 

n=l  m=l 

Now  consider  an  environment  with  a  direct  path,  n=l,  with  arrival  angle  near 
the  horizontal,  a1>am  for  m^l.  and  *£>0.  The  arrival  of  the  direct  path  near  the 
horizontal  allows  the  bearing  to  be  estimated  from  as 

'?=sin“1[|r1j  (A-27) 

Because  aj>am>  *11  m^l,  the  n=m=l  term  of  (A-26)  is  the  largest,  given  by 
2 

a^gfr-Tj).  Note  that  this  is  exactly  the  peak  that  would  occur  if  there  were  only 
a  direct  path. 

The  effect  of  the  other  terms  of  (A-26)  will  depend  upon  their  proximity  to  the 

peak  at  Tj.  If  •yQ— >^+1^^,  then  the  term  will  tend  to  modify  the  peak  atr=i^, 

biasing  it  in  the  direction  ofr=y -y  +t  and  broadening  it.  Since  there  is  no 

n  'm  n 

reason  for  the  terms  of  (A-26)  near  tst^  to  be  symmetrical  about  then  the 
multipath  structure  may  bias  the  peak  away  fromr=T1.  The  broadening  of  the 
peak  can  be  interpreted  as  a  reduction  in  the  effective  bandwidth  of  the  signal 
process.  From[l],  Appendix  IV,  the  standard  deviation  of  the  bearing  estimate 
is  inversely  proportional  to  bandwidth  to  the  3/2  power,  while  inversely  propor¬ 
tional  to  signal-to-noise  ratio.  Thus,  a  given  percentage  decrease  in  the  effec¬ 
tive  bandwidth  due  to  multipath  will  degrade  tracker  performance  more  than  an 
equivalent  decrease  in  SNR. 

Other  terms  of  (A-26)  that  are  peaked  at  values  of  mot  near  will  produce 
additional  signal  related  peaks  in  the  mean  weight  vector.  These  peaks  will,  in 
turn,  be  biased  and  broadened  by  other  multipath  arrivals  with  a  peak  in  the 
same  vicinity.  It  would  be  possible,  given  an  unfortunate  multipath  structure, 


for  one  of  these  other  peaks  to  be  larger  In  amplitude  due  to  the  combination  of 
several  paths.  However,  by  use  of  a  windowing  scheme  as  described  in 
Section  2.3  for  use  with  multiple  targets,  it  will  be  possible  to  track  on  the 
peak  atr=r1  provided  the  window  is  narrow  enough  with  respect  to  the  separation 
of  the  various  peaks. 

It  is  difficult  to  further  generalize  the  effects  of  multipath  on  the  adaptive 
tracker  structure,  since  the  effects  depend  very  much  on  the  particular  multipath 
structure  present.  Potential  effects  include  biasing  of  the  bearing  estimate. 
Increased  estimation  variance,  and  ambiguity  in  the  weight  vector  as  to  the 
correct  peak  for  tracking.  However,  these  factors  affect  the  adaptive  tracker 
in  very  much  the  same  way  they  do  conventional  trackers,  and  their  sensitivities 
to  multipath  will  be  similar. 
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THE  EFFECTS  OF  ARRAYS  ON  THE  ADAPTIVE 
TRACKER  WITH  STATIC  TARGETS 

[i  o] 

Throughout  the  first  two  phases  of  the  adaptive  tracker  study  ,  the 
inputs  to  the  adaptive  tracker  have  been  modeled  as: 


x(t)  =  s(t)  +  nx(t) 

(B-l 

d(t)  =  s(t-r)  +  nd(t) 

(B-2 

where  s(t)  is  the  plane  wave  signal,  n^t)  and  n^(t)  are  random  noise,  uncorre¬ 
lated  with  each  other,  and  r  is  the  delay  in  the  signal  due  to  the  offset  of  the  tar¬ 
get  from  the  steering  angle.  If  the  sensors  consist  of  two  omnidirectional  hydro¬ 
phones,  each  driving  one  of  the  tracker  inputs,  then 

r  =  fsin  9  (B-3 

c 

where  d  =  distance  between  sensors 
c  =  speed  of  sound 

0  =  angle  between  plane  wavefront  and  line  passing 
through  the  hydrophones 

This  model  may  also  be  applied  to  the  more  general  situation  of  Figure  Bl, 
in  which  a  half  array,  steered  in  the  general  vicinity  of  the  target,  drives  each 
tracker  input.  In  this  case,  d  must  be  replaced  by  the  distance  between  the  phase 
centers  of  the  half  arrays.  However,  as  the  mis-steering  angle  deviates  from 
zero  (that  is,  the  target  bearing  differs  from  the  steering  angle),  the  frequency 
response  of  the  half  array /be  amforme  rs  change,  modifying  the  signal  spectrum. 
The  model  given  in  (B-l),  (B-2),  and  (B-3)  has  therefore  been  used  as  a  good 


approximation  when  the  mis -steering  angle  is  small.  The  analysis  that  follows 
will  show  that  the  model  used  above  applies  for  a  very  general  class  of  arrays 
even  when  the  mis-steering  angle  is  not  small,  provided  the  spectrum  of  s(t) 
is  appropriately  modified. 

Let  each  half  array  consist  of  K  hydrophones  arbitrarily  arranged  in  space 

th 

as  shown  in  Figure  B2,  where_Ym  is  a  vector  from  the  origin  to  the  m  hydro¬ 
phone.  If  JJt  is  a  unit  vector  in  the  direction  of  the  target,  andJ2g  is  the  steering 
angle  of  the  array,  then  the  output  of  the  hydrophone  can  be  written  as 

Zm<*>  *  »  ['  -4-VX  - +  °i  [*  -&A]  (B-4) 

where  the  n^(t)  are  uncorrelated  with  each  other.  Now,  let  m-1, 2, ...  K  be  the 
left  half  array  elements  with  shading  coefficients  am,  and  m*K-4, ...  2K  be  the 
right  half  elements.  Then  the  outputs  of  the  two  half  arrays  are 

*<*>  ■  »i  <B-5> 

2K  ) 

d(t) "  rn5n  ai  ^  +  n  '‘"eVVj  (B'6) 

The  mean  weight  vector  of  the  M  tap  adaptive  filter  is  the  discrete  Wiener 
filter,  given  by 

E  [2L  (s.s)]  (B-7) 

where  R^  is  an  MxM  matrix  with  the  p,  q  element  given  by 

'Vp  v E  x  (,"pv  <B-8' 

where  T  Is  the  adaptive  algorithm  sample  rate, 
s 


and  Is  a  K  dimensional  column  vector  with  the  q**1  element  given  by 

(Rdx)q  3  E[x(t-qTg)  d*(t»  (B-9) 

Now,  let  p (t)  and  P  (t)  be  the  normalized  covariance  functions  of  the  signal 

3  n 

and  noise  (assumed  wide  sense  stationary) , 

PJt)  3  4  E  fx(t)  x*(t+r)J  (B-10) 

x  rx 

with  Px  the  power  In  x(t).  Denote  the  power  spectral  density  of  the  signal  and  noise 
respectively  as  S  («)  and  S(cj).  Then,  using  (B-5)  and  (B-6), 

3  u 

K  2K  *  r  ,  .1 

<Rdx>q  =  PS  £  k§+1  %  Ps  [qTS  *f*k>  ^“^>1  <B-n> 


Writing  this  in  terms  of  the  power  spectral  density  of  the  signal,  where 


00 

ps 

—  00 


(B-12) 


gives 


(B-13) 


du> 


B-4 


Similarly, 


^Rxx\>i<J  ~  2* 


fs  (1^mph,)ts  f£  f  ..-igsJu-iy 

-<C  8  1=1  1  k=K+l  K 


27rXSn(w)eja,<P"q)Ts  2  |ai|  l  dw 

-00  1=1  [  1  J 


du 


(B-14) 


Equations  <17)  and  (18)  can  be  simplified  by  defining 


K 


HdM  -  E  (VJV 

Q  k=l  K 


(B-15) 


and 


2K 


H  (u>)  =  g  a.  e^k  <V*t> 

x  k=K+l  K 


so  that 


(B-16) 


00 

•*00 


(B-17) 


and 


1  |  80  K 

<‘Wp,q  E 

-•  I  Zoo  1=i 


at  dw 


(B-18) 


Note  that  the  model  of  (B-l),  (B-2),  and  (B-3)  with  a  single  omnidirectional 
hydrophone,  K=l,  would  give 


HXM  =1 


The  model  can  then  be  interpreted  aa  replacing  the  half  arrays  and 
beamformers  with  the  linear  filters ,  Hx(w)  and  H^(cj),  as  shown  in  Figure  B2. 
Suppose  that  the  array  geometry  is  constrained  so  that 

— k+M  *  — k  +  — o  (I 


and  that  the  shading  coefficients  satisfy 


=  \  k  ®  1,2, ...  ,K 


(B-22 


Figure  B2.  Equivalent  Model  for  Split  Array  Tracker 
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Then  (B-16)  cm  be  written 


H  M  = 
x  ' 


E  ^kX--“«> 

k=K+l  K 

E  JAX  -b,> 

K 


or 

Hx(w)  =  Hdfu)^o  (iJs"^t) 
In  this  case, 


H  (w) 
x' 


A,  +. 


2  ^  (Ug-Ut)dw 


and 


'Vp,q '  s/V'  M2  ‘1"tTir A<“>  E, 

-ao  -oo  1=1 


(B-23) 


(B-24) 


(B-25) 


This  is  exactly  the  same  form  as  the  case  with  two  omnidirectional  hydrophones 
if  the  signal  spectrum  is  replaced  by  Sg(w)  |iyw)|2  and  the  noise  power  scaled 


by 

K  ,  i  2 


In  this  equivalent  model,  the  delay  between  the  outputs  of  the  two  hydrophones  is 


T  a 


COS  4> 


(B-26 ) 


B 


where  4>  is  the  angle  between  2L  and  (J^-J^).  This  equivalent  model  is  shown 
in  Figure  B3.  Thus,  the  analyses  of  [1]  and  [2]  can  be  applied  when  the  array 
meets  B-21)  and  (B-22). 

A  number  of  array  configurations  of  interest  meet  these  constraints.  For 

example,  let  the  full  array  consist  of  uniformly  spaced,  unshaded  elements  on  a 

straight  line,  with  spacing  d  .  Then 

s 


sin(^ud) 

H,(cj)  = - f - e  '2  M 

*  sin(jwo$ 


(B-27) 


and 

Md 

T  — - -sin  0  (B-28) 

c 

with  the  angle  between  the  array  and  the  plane  wavefront.  Incorporation  of  non- 
uniform  shading  satisfying  (B-22)  will  change  Hd(o;),  but  notr.  More  generally, 
any  two  identical  arrays  with  identical  shading  and  arbitrary  relative  position 
satisfy  the  constraints  in  (B-21)  and  (B-22),  and  can  be  modeled  using  the 
omnidirectional  hydrophones . 


B-8 


"d<‘> 


Figure  B-  3.  Equivalent  Model  for  Split  Array 
Tracker  in  Special  Case 
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APPENDIX  C 


EFFECTS  OF  ARRAYS  ON  DYNAMIC  TRACKING 

In  the  second  phase  of  the  Adaptive  Tracker  Study ,  [  2]  the  performance  of 
the  adaptive  tracker  in  the  presence  of  dynamic  targets  was  investigated  when  the 
tracker  inputs  were  provided  by  two  omnidirectional  hydrophones.  This  was 
done  under  the  assumption  that  the  array  effects  could  be  introduced  later  by 
modifying  the  input  spectra.  Appendix  B  of  this  report  showed  that  his 
assumption  was  true  for  static  targets.  This  section  considers  the  effects  of 
split  arrays  in  the  dynamic  case. 

As  in  the  omnidirectional  hydrophone  case,  when  the  target  is  dynamic,  the 
weight  vector  is  continually  changing  with  the  peak  of  the  time  domain  weight 
vector  tracking  the  true  delay.  It  would  be  expected  that  eventually,  the  peak 
would  lag  a  linearly  varying  delay  by  a  fixed  amount ,  corresponding  to  the 
steady  static  lag  of  a  first  order  system.  This  is  complicated,  however,  by 
the  fact  that  a  linearly  moving  target  will  eventually  leave  the  mainlobe  of  a 
beam ,  either  requiring  that  the  beam  be  steered  or  the  tracking  operation 
handed  over  to  another  beam  ( for  a  discussion  of  this  process ,  see  Section 
2.4).  In  either  case,  this  introduces  a  discontinuity  in  the  position  of  the 
target  relative  to  the  maximum  response  axis  (MRA)  of  the  tracking  beam. 
Further ,  if  the  weights  are  reinitialized  to  reflect  the  bearing  estimate  already 
computed,  this  introduces  a  discontinuity  in  the  weight  vector. 

The  goal  of  this  analysis  is  to  determine  the  main  steady  state  weight  vector 
of  the  adaptive  tracker  with  split  arrays  and  a  broadband  target  with  linearly 
varying  time  delay.  However,  since  a  target  with  linearly  charging  delay  will 
eventually  leave  the  mainlobe,  the  steady  state  mean  weight  vector  for  this 


C- 


t 


target  is  zero.  This  difficulty  was  avoided  by  Introducing  a  periodicity  into  the 
time  delay  input  as  shown  in  Figure  C-l.  This  periodic  waveform  allows  the 
target  to  remain  in  the  beam  so  that  a  steady  state  mean  vector  can  be  obtained 
at  the  point  when  the  target  is  about  to  leave  the  mainlobe,  as  shown  in 
Figure  C-l.  The  periodicity  is  a  mathematical  artifice  only,  and  does  not 
represent  any  aspect  of  the  physical  problem.  This  periodicity  is  shown  to 
produce  a  separate  term  in  the  expression  for  the  lag  in  the  time  delay  estimate, 
so  that  its  effect  can  be  removed  by  neglecting  that  term. 

From  reference  [2] ,  the  mean  weights  of  the  adaptive  filter  are  given 


EH  ■  *  £  f-“Rxx<n-'l-1)]D"k'lRdx'k>  <0-W 


where 


W(n)  =  n  time  domain  weight 
p  =  feedback  coefficient 

Rxx(n)  *  covariance  matrix  of  adaptive  filter  data  vector 
th 

on  n  iteration 

Rgx(k)  *  cross  correlation  between  data  vector  and  desired 
input 


Consider  the  outputs  of  two  split  array  beamformers  operating  on  a  uniformly 

spaced  line  array  with  a  dynamic  target.  The  signal  propagates  across  die 

array  with  a  time  varying  delay,  a(t) ,  between  adjacent  hydrophones.  The 

beamformers  apply  a  delay  r  between  adjacent  hydrophones  and  a  bulk  delay, 

s 


C-2 


Therefore,  the 


th 

q  element  of  Rdx(n)  is 


r  t  M  M  2 

E  Jx(t)d(t+q6)j  *  52  L  <rs 


j-q6+  |m-hJ  fa(t-Hj6)-TsJ  - 


In  order  to  use  the  triangular  input  given  in  Figure  Cl,  let 


(C-8) 


p(t)  *  a(t)-ra  =r(t-ST)-rp 


where 


<k  <(S+l>g-  and  r  = 


(C-9) 


Here,  T  is  the  period  of  the  sawtooth,  2 rp  is  the  peak  to  peak  amplitude  of 
the  sawtooth,  and  S  is  the  index  of  the  cycle  of  the  sawtooth  (i.e. ,  S=2 
indicates  the  second  cycle) . 

Now,  it  is  assumed  that  the  signal  to  noise  ratio  is  sufficiently  low  that 


R^Ct)  -  E  [x(t)xT(t)  J  *  Mon2  I  (0-10) 

2 

where  M <7  is  the  noise  power  after  beamforming.  Then  to  evaluate  (C-l), 
n 

(C-8)  and  (C-10)  are  used  with  the  substitution  t*k6 ,  so  the  mean  weight  is 

eKH  n’k’1  J-  g  *.%,  }-o*+H[p(k«^] 


-Jp(k8) 


(C-ll) 


Now 


ST. 


p(k8-tq6)  ■  rSflc+q^y)  -  Mrp  *  p(k5)  +  r<Sq 
so  letting  p(k5)  *  p(k). 


(C-12) 


M  M  M 


[»  .  A  m  *  ' 

W  (n)l  -pff2£  £  £  (1-hM<t  2) 

•  q  J  3  i=l  j=l  k=0  n 


2  n-k-i  _  | 


pa  j(M+i-j)p(k) 


(M+i)rqS  -q5 


(C-13) 


It  Is  now  assumed  that  the  signal  has  an  exponential  correlation  function. 


PS(T)  =P'TI  ,  P€  (0,1) 


(C-14) 


so 


L  q  J  s  i=l  j=l  k=0  n 


(C-15) 


Now  suppose  that  (M-H)r«  1,  so  that  (M+i)rq  can  be  neglected  with  respect  with 
to  q.  Therefore,  (C-15)  can  be  approximated  as 

P 1  6  1  (c-i6) 

rT 

Further,  using  r=2i^/T  requires  xp  =  so  that 
p(k)  -  r6  [k-(S+|)j] 


(C-17) 


Using  this,  It  is  possible  to  evaluate  the  inner  summation  in  (C-16)  by  taking 
advantage  of  the  periodic  nature  of  p(k).  Assume  that  the  period  of  the  input  is 
an  integer  number  of  iterations,  i.e. ,  T/S  is  an  integer,  so  that  (C-16)  can 
be  written 


M 


0  M  M 

»°s  £  £ 

3  i=l  j«l 

p5|(M-H-j)r[k-(e^f] 


3=0  k*sT  n 

—<5 


(C-18) 


.1,T 


Letting  v=k-(S+J)|-  in  (C-18)  gives 


r  ~\  2  &  &  |W>T  ,n-l-(S^)£  T/26  , 

[wqw].„<r,2  E  r  [g  a-vMcJb 


V-_  26 


s|(M-*i-J)  rv-q 


(C-19) 


These  two  sums  can  be  calculated  separately. 


(WM<7Q2)n  1  26 


2  -(n-D-y 


T 


(C-20) 


As  n— »  -  ,  the  first  sum  goes  to 


lim 

00 


(n-1 4 

£ 

3=0 


2  D-1'(9+^ 
(1-hM  aa£)  26- 


2  -T/26 


«  T/d 


(C-21) 


This  assumes  that  1 1-uMa  2 / 

atl<WSNR-  that  this  Lm  is  ^  ^  ^  °^tion  fo*  stability 
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(C-23) 
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1  .  p(M+H)  ir 


1-p.Ma  2 
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(C-24) 


This  must  now  be  summed  on  i  and  j  as  in  equation  (C-19).  To  do  this,  let  i=j 
in  the  denominator,  so 


p  (M+i-j)  5r  xpU6r 


which  is  approximately  true  if  Mbr<5  <<  1.  Then  (C-24)  reduces  to 


2,T/25  (X-^mO 
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+  1  _Mr5  2"1 

1-P  (1-tiMo^  ) 


1-p  Mr5  (1-^M  an2f 1  (C-25) 


Now,  approximating  the  sum,  (C-19),  by  letting  i*j  in  the  second  and  third  terms 
of  (C-25)  gives 

E  [Wq(n)J  =  h-  <rs2  [K3L'a^Mo-n2)’q/Mr  -  K yq  -  Kg'p-^J  (C-26) 


where 
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(C-28) 


The  time  delay  estimate  will  correspond  to  the  peak  of  the  mean  weight 
function,  so  assume  that  q  is  a  continuous  parameter  as  in  Appendix  G  of  [2]. 
In  order  to  obtain  die  peak  relative  to  the  true  delay  (i.e. ,  the  lag)  let 


q  =  -  p  ,  p  St  0 


Then  (C-26)  can  be  written 


(C-29) 


Mr 

E  [Wq(n)l  =  hl«ts2  [^(1-hM^2)  -KgP  “6P  -  K3P6p]  (C-3 
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'26  . 


Kx  = 


K  =  PMrT/2K' 

2  * 

K3  =  P"MrT/2K’3 


(C-31) 


Noting  that 


-2.  2.  „  2,  flMr«l 

g  =*  P  (l-nM<rn  )  |^(l-v.M<rn  )  P  J 


(C-32) 


C-10 


the  mean  weight  can  be  further  approximated  as 


[vn)]  *  >-s2[ki 


(nxM<rnA) 


p/Mr 
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(C-33) 


To  determine  the  peak  of  this  function,  differentiate  with  respect  to  p  and 
set  to  zero,  yielding 


p  =  Mr _ KHn(l-pM^., 

£n  (1-jaM<t2)  -£npMr<S 


(C-34) 


as  the  lag  in  taps  in  the  peak  of  the  mean  weights  in  comparison  to  the  true  delay. 
2 

If  (j.M<rQ  «  1  and  for  an  exponential  correlation, 


with  pMr<s  «  i-Mrb<S,  then  the  lag  reduces  to 


Lag- 


Mr 


-  bMc6 


|  .69  +2  2n 


M(l-e  br^  (l.+T/26))' 
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(C-35) 


This  can  be  compared  to  the  analogous  results  without  the  array  present  given 
in  Appendix  G  of  [2]  it  it  is  noted  that  Mr,  the  rate  between  array  centers, 
corresponds  to  c  in  the  reference.  The  first  term  is  identical  to  the  result 
obtained  without  the  use  of  arrays,  that  is,  with  a  pair  of  omnidirectional  sensors. 
The  second  term  represents  the  additional  lag  introduced  by  the  array.  It 
appears  that  the  third  term  Is  a  reduction  in  the  lag  due  to  the  periodicity  of  the 


input,  anH  is  therefore  a  product  of  the  mathematical  artifice  introduced  to  solve 
the  problem,  not  of  the  device  itself.  The  lag  in  taps  can  therefore  be  approxi¬ 
mated  as 


Mr 

Lag  *  o 

liMo-n  -  bMrfi 


.69+2  in 
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Note  that  this  reduces  to  previous  results  for  M=l.  This  can  be  written 


Mr 


-bMrS 


J  .69  +  2in  M  +  2£n 
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Lag 


Now,  from  (C-9), 


(C-38) 


where  2tp  is  the  peak  to  peak  amplitude  of  the  sawtooth.  Taking  tp  to  be 
the  delay  associated  with  the  target  at  the  edge  of  the  mainlobe,  i.e.,  when 
the  deviation  from  the  MRA  equals  the  beamwidth,  assuming  a  distance  Md 
between  phase  centers,  gives 


2t 

t  =  — ^  sin  B 
p  r  w 


(C-39) 


with  B  the  beamwidth.  Then 
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2Md 

rc 


sin  B 
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(C-40) 


and 


br5(l  +  T/25)  =  b(r«  +  —  sin  B  )  (C-41) 

c  w 

Note  that  r<5  is  the  change  in  delay  of  the  target  in  one  algorithm  iteration 

while  d/c  sin  B  is  the  time  delay  associated  with  a  target  at  0  =  B  .  In 
w  w 

practical  applications  r6  <<  d/c  sin  B  ,  so 
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The  third  term  in  the  bracket  will  be  neglible  in  most  cases,  which  can  be 
shown  as  follows.  Assume  that  the  array  design  frequency  (the  frequency  at 
which  the  hydrophones  are  one-half  wavelength  apart)  is  b,  so  that  the  array 
is  designed  for  the  highest  frequency  in  the  band.  Then 


and 

bT  8inBw  =  fsinBw  (C"43) 

if  the  beamwidth  is  less  than  30° ,  then  sin  B  «  B  ,  and 

w  w 

M  mb 

"  *>"  ®.  -  -T1  (c-44) 


Clearly  M  and  Bw  are  inversely  proportional.  A  reasonable  approximation, 
based  upon  the  beamwidth  of  a  line  array,  is  that 

Bw  -  m  (C-4S> 

with  L  the  array  length.  The  half  array  beamwidth  chosen  above  is 
(M-l)  c/2b,  where  \  =  c/b.  Then 


MBw  _  M  50  M 
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The  third  term  in  the  bracket  is  then 
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which  is  virtually  unity  for  all  values  of  M  >1.  Therefore,  the  lag  is  taps  is 
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In  order  to  interpret  this  as  a  bearing  lag,  assure  that  the  bearing  angle  is 
near  enough  to  broadside  so  that  the  bearing  9,  is  related  to  delay,  r,  by 


where  c  is  the  speed  of  sound  and  D  the  distance  between  array  phase  centers. 
For  the  array  considered  here, 


D  =  Md 


and  the  lag  In  seconds  is  (Lag)  <5  .  Therefore,  the  lag  in  bearing  is 
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APPENDIX  D 


ESTIMATOR  STRUCTURES  FOR  SIGNALS  WITH  BOTH  BROADBAND  AND 
NARROWBAND  COMPONENTS 

As  pointed  out  in  the  introduction,  one  of  the  motivations  for  use  of  an 
adaptive  tracker  is  that  is  has  the  potential  to  track  signals  using  both  broadband 
and  narrowband  energy  simultaneously.  Reference  [1]  developed  estimator 
structures  for  use  when  either  broadband  or  narrowband  energy  only  were 
present  in  the  radiated  noise  of  the  target.  References  [1]  and  [2]  then  analyzed 
the  performance  of  these  two  estimators  in  the  presence  of  targets  that  are 
broadband,  narrowband,  or  a  combination  of  the  two.  However,  a  means  of 
combining  the  two  estimates ,  or  performing  a  different  estimate  when  the 
target  has  both  broadband  and  narrowband  components  was  not  considered. 

In  this  section,  estimator  structures  that  are  tailored  to  targets  with  both 
broadband  and  narrowband  energy  are  developed .  Two  of  the  structures 
are  based  upon  a  maximum  likelihood  estimate  of  the  target  bearing  given 
certain  assumptions  on  the  weight  vector  statistics.  However,  several  of 
these  assumptions  have  not  been  proven  rigorously,  and  are  made  primarily 
for  analytical  tractability .  Therefore  the  resulting  structure  cannot  right¬ 
fully  be  called  the  maximum  likelihood  estimator.  The  estimator  does  combine 
the  elements  of  the  broadband  and  narrowband  trackers  described  in 
reference  [  1] . 

In  the  adaptive  tracker ,  the  bearing  estimate  is  extracted  from  the 
weight  vector.  Assume  that  the  weight  vector  is  a  complex,  gaussian  random 
vector , 

W_  =  [W(0),  W(l),  ...W(M-1)]T  CD-I) 


D-l 


with  mean 


r  -l  -jw  Oct  -  -  sin  0)  . 

M0O  a  E  [W^)]  *  ae  °  80  +  bpCkr  --sine)  (I>2) 

w  sc 

The  function,  p(t),  is  an  arbitrarily  normalized  function,  corresponding  to 
the  shape  of  the  weights  due  to  the  broadband  component.  Also, 

u>0  =  radian  frequency  of  narrowband  component 

T  3  algorithm  sample  rate 
s 

d  =  distance  between  split  array  phase  centers 
c  3  speed  of  sound 
9  =  bearing  of  target. 

The  gaussian  property  of  the  weights  can  be  argued  based  upon  the  large  number 
of  iterations  of  the  weight  update  algorithm,  while  the  structure  of  the  mean 
weight  is  from  references  [1]  and  [2].  Now,  assume  that  the  weight  vector 
has  arbitary  covariance  matrix,  R  ,  which  is  Independent  of  target  bearing,  hi 

W 

[2] ,  the  variance  of  the  frequency  domain  weights  are  shown  to  be  independent 
of  target  bearing,  so  it  seems  reasonable  that  this  property  should  hold  for  R^. 

Given  these  assumptions,  the  log-likelihood  function  of  the  weight  vector  can 
be  written  as 

*(W|9, a,b)  =  Cx  +  1/2  (tf  -  Slw)+  R"1  (W  -  ®w)  (D-3) 

where  is  a  constant  independent  of  6  and  the  superscript  +  denotes  conjugate- 
transponse.  The  maximum  likelihood  estimate  of  6  is  that  value  of  6  which 
maximizes  the  second  term  of  (D-3).  Expanding  (D-3)  gives 

*<W|9,a.b)  3  Cx  +  1/2  WV1  W  + 1/2  m^  R^1  -  Re  [w+R"^  W) 
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Hie  first  term  of  (D-4)  does  not  depend  upon  the  variables  to  be  estimated,  so 
it  need  not  be  considered. 


In  order  to  proceed  further,  it  is  assumed  that  the  weights  have  equal  variance 
and  are  uncorrelated,  that  is, 

R  *  <t3.  I 
w  W 

The  assumption  of  equal  variance  has  been  demonstrated,  at  least  for  narrow- 
band  plus  white  broadband  signals  in  [l]  and  [2] .  However,  the  weights  are 
probably  not  uncorrelated  when  the  sample  rate  is  several  times  the  Nyquist 
rate,  as  will  be  the  case  in  practice.  However,  the  assumption  does  yield  a 
physically  satisfying  result,  even  though  it  is  not  truly  the  maximum  likelihood 
estimator. 

With  the  assumption,  the  second  term  of  (0-4)  can  be  written 


a;lv  -  _i  2  if  { aW, v,-j  «<■>•>  <»-»> 

2g  K"0  i 
w 

+  2ab  p(kTs-^sin0)  cos  Ju>o(kTs-£  sin0) 

If  the  length  of  the  adaptive  filter  is  long  in  comparison  to  the  period  of  the 

sinusoid  and  the  correlation  time  of  the  narrowband  component  and  the  delay  is 

not  near  either  end  of  the  tipped  delay  line  (that  is  d/c  sin  0  is  not  close  to  zero 

or  (M-l)T  ),  this  is  essentially  independent  of  8.  The  third  term  of  (D-4)  is 
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Re 
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£  W(k) 
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+  bp(kTs--  sin0) 


-jus  (kTs-§  sin8) 
ae 


(0-6) 


Then  the  joint  estimate  of  9,  a,  and  b  is  that  ( 9,a,b)  which  minimizes 


1  M_1  / 

J(6,a,b)  » —y  j  -  [a2 +b2  p2(kTa -^sinS) 


+  2abp  (kTs  -  ~  sln8 )  cos 
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+  Re 
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Now,  let  W^(m)  be  the  FFT  of  the  weight  rector. 


M-l  ,2ir  K 

W-(m)  =  £  W(k)  J  M (D-8) 
1  k=0 

and  let  W^(max)  be  the  frequency  domain  weight  with  the  largest  magnitude 
(containing  the  signal).  If  the  FFT  resolution  is  sufficiently  fine  that  (he  center 
frequency  of  the  largest  bin  is  a  good  approximation  to  wQ,  than  (D-7)  can  be  well 
approximated  as 


J(&,  a,b)  *  — iy  -Ma2  -  Ct  b2  -  2C2ab  (I>9) 

2a  ' 
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d  M-l  d 
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where 
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(D-10) 


Under  the  assumptions  made  following  equation  (D-5),  C,(  and  C2  are  essentially 
independent  of  6. 

Now, 


9  J 
9a 


j  j“2Ma  -  2C2b  +  2 1 W  (max)|  008(4^-  ^  sin  9) 
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E  W(k)p<lcTs-*o-sin0) 
k=o  ) 


(D-ll) 


(D-12) 


Setting  these  to  zero  yields  the  maximum  likelihood  estimate  of  a  and  b  provided 
9  also  minimizes  (D-  9).  Then 
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and 
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(D-14) 


The  minimization  over  9  cannot  be  carried  out  explicitly,  and  must  be  done 
numerically  using  (D-13),  (D-14),  and  (D-9).  This  could  be  done  by  the  follow¬ 
ing  procedure: 

(a)  Select  an  initial  bearing  estimate,  e0 
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(b)  Extract  the  magnitude  of  the  largest  FFT  bin,  |w^  (max)  | ,  the 
phase, and  its  frequency,  Wq. 

(c)  Numerically  minimize 


and 
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(D-16) 


-  j|  [  Wf(max)|  coe^-w^sinSj)  ( 


(D-17) 


Here  C  and  C„  are  constants ,  essentially  independent  of  81  given  by  (D-10),  and 
9;  is  the  value  of  the  estimate  on  the  Ith  iteration  of  the  minimization  procedure. 
Note  that  the  data  extracted  from  the  weights  remains  the  same  during  the 


iterations  on  i. 


This  structure  can  be  simplified  by  noting  that  by  substituting  (D- 16)  and 


(D-17)  in  (D-15)  J  can  be  written  in  the  form 

•J  (6j)=  Q^Qj)  +  +  2Q3f1(9i)f2(0.) 


(D-18) 


where 

f  1  (©i)  =  |  Wj(max)  |  cos(<()w-  sinQj) 


(D-19) 
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and 
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Q3  =  C“  (M+l)  -  C^C^l) 


(D-21) 

(D-22) 

(D-23) 


so  that  Q1?  Qg,  and  Qg  are  constants  essentially  independent  of  9^  Then  step  (c) 
in  the  numerical  procedure  becomes 


(cl)  Numerically  minimize 

J$i)  =  Q^j)  +  Q2f^  (@j)  +  2Q3f1(@i)f2(@i)  (D-24) 

where  is  the  value  of  the  estimate  on  the  ith  iteration  of  the  numerical  proce¬ 
dure.  This  leads  to  the  simpler  structure  shown  in  Figure  Dl. 

An  alternate  approach  is  possible  if  the  broadband  sign al-to- noise  ratio 
(BSNR)  and  narrowband  signal  to  noise  ratio  (NSNR)  are  known  a  priori  or  are 
available  from  the  sonar. 
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Figure  Dl.  Combined  Adaptive  Tracking  System  with  Unknown  SNRs 
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From  [  2] ,  the  kth  element  of  frequency  domain  weight  vector  in  steady  state 
for  a  signal  consisting  of  a  white  broadband  source  with  a  single  line  component 
is,  neglecting  FFT  windowing  effects 


Tg  =  sample  interval 

A  =  signal  delay  between  phase  centers  in  sample  intervals 
(assumed  to  be  an  integer) 


(D-25) 


(D-26) 
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Now,  if  the  narrowband  signal  is  bin  centered,  or  approximately  so,  then 

|  G(k,a)o)  j  ^  58  M2j  ^  (D-27) 

where  5^  is  the  Kronecker  delta  and  J  is  the  index  of  the  bin  containing  the 

narrowband  component.  If  windowing  effects  are  negligible,  as  they  will  be  in  a 

well-designed  tracker,  then  the  time  domain  weight  vector  can  be  obtained  by 

til 

inverse  transforming  (D-25).  If  (D-27)  is  valid,  then  the  a  time  domain 
weight  is  approximately  given  by 
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This  result  Indicates  that  the  mean  weight  used  in  the  derivation  of  the  combined 
estimator  is 


Mw(k)  = 
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When  these  coefficients  replace  a  and  b  in  the  derivations,  the  combined 
estimator  must  minimize  the  functional 
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W^iere  Y  NB  =  °NB  YBB  =  °BB  ^  °N  ’  Hence»  S^n  a  priori  or  external 

knowledge  of  vNB  and  YBfi,  step  (c)  in  the  minimization  procedure  becomes 


(c2)  Numerically  minimize 


Jl(6i)  (MY^  +  l)  <rBB  +  l)  I  wf  (max)|  008  (*w  ““of  sin  6i) 
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(D-31) 


k=0 


This  estimator  structure  is  shown  schematically  in  Figure  D2. 
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This  reduces  to  the  broadband  or  narrowband  estimator  when  the  signal  is 

o 

purely  narrowband  or  purely  broadband.  When  <r_,_  =0,  then  =  0  and 

DD  Bfi 

y 

Jl(®)  =  1  lWf  (max)l  cos  (*w  8in  (D-32) 

NB 

which  is  minimized  by  choosing 

A 

6  =  arc  sin 

which  is  exactly  the  procedure  used  for  narrowband  signals.  If  =  0,  then 

NB 

the  procedure  reduces  to  the  minimization  of 
M-l 

J1(0)  =  J  w(k>p<kT3  *  f  sin$)  (D-33) 

k=l 

This  is  just  a  minimization  of  the  interpolated  time  domain  weights  using  the 
interpolation  function  p(t).  This  is  the  broadband  estimator  if  p(T)  is  taJken  tb 
be  the  sin(x)/x  function.  This  result  is  a  strong  justification  of  the 
narrowband  and  broadband  estimators  originally  described  in  [1] . 

One  drawback  of  these  estimator  structures  is  that  they  are  highly 
parametric.  They  not  only  require  knowledge  of  (or  estimation  of)  the  signal 
powers,  but  assume  a  priori  that  the  signal  consists  of  a  broadband  component 
and  a  single  line  component  in  a  non- directional  background.  They  may, 
therefore,  be  very  sensitive  to  deviations  of  the  actual  environment  from  the 
assumptions,  i.e.,  multiple  target  lines,  interferences,  etc. 


«>W 
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A  completely  different  philosophy  of  providing  estimates  for  targets  with 
broadband  and  narrowband  components  is  to  combine  the  outputs  of  the  broadband 
and  narrowband  estimates  according  to  some  weighting  scheme.  The  weights 
are  chosen  according  to  a  criterion,  such  as  minimum  variance  of  the  combined 
estimate.  Suppose  initially  that  the  two  estimates,  and  §Bfi,  were  uncor¬ 
related  estimates  of  the  bearing  of  a  single  target  from  the  narrowband  and 
broadband  estimators ,  respectively.  It  can  be  shown  that  the  linear  combina¬ 
tion,  6„,  of  the  two  that  yields  the  same  mean  and  minimum  variance  among  all 

v 

linear  combinations  is  given  by 
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where 
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=  variance  of  8xrQ 
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If  0.  is  the  variance  of  the  combined  estimate,  then  it  can  be  shown  that 
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The  approach  is  easily  extended  to  combining  more  than  two  estimates.  There¬ 
fore,  if  the  target  had  multiple  lines,  a  narrowband  estimate  could  be  extracted 
from  each  and  combined  in  this  way. 
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It  has  been  assumed  that  the  estimates  come  from  a  single  target,  while,  in 
fact,  multiple  targets  and  interferences  may  be  present.  A  reasonable  assump¬ 
tion  is  that  if  two  estimates,  or  the  mean  of  two  estimates,  are  close  together, 
then  they  come  from  the  same  target  and  can  be  combined  in  this  way.  The 
combiner  may  therefore  be  implemented  (for  a  single  broadband  and  a  single 
narrowband  estimate)  as  follows. 

a.  Compute  a  running  sample  mean,  M '  and  M  *  ,  and  a  running 

°NB  eBB 

sample  variance,  V:  and  V'  ,  for  the  narrowband  and  broad- 
A  ®NB  A  SBB 
band  estimates,  e^g  and  eQB 

b.  If  jM£  -  Ma  I  <  g,  where  £  is  a  threshold,  then  combine 
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and  eBB  to  get  ec  as  follows: 
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It  should  be  clear  that  the  sample  variances  in  (D-36)  are  random  variables,  and 
that  their  fluctuations  can  increase  the  variance  of  e  to  the  point  where  (D-35) 

C 

is  no  longer  true. 

To  this  point,  it  has  been  assumed  that  and  §BB  are  uncorrelated. 
However,  since  they  are  extracted  from  the  same  data,  the  weight  vector,  this 
is  unlikely.  Of  course,  it  is  still  possible  to  use  the  above  combiner  because  of 
its  simplicity,  but  the  variance  of  §c  will  be 
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8NB  ^BB  =  E  [  (<r®NB  "  E(6NB))  (®BB  '  E(0BB))  * 

That  is,  p  is  the  correlation  coefficient  between  and  §fiB.  Note  that  increas¬ 
ing  correlation  always  increases  the  variance  of  §  .  This  can  be  written 
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Consequently,  if  either. (D-40)  or  (D-41)  are  true,  the  combiner  cannot  reduce  the 

variance  beyond  that  of  the  individual  estimates.  Figure  D3  shows  T^  and  Tg 

plotted  as  a  function  of  (<r£  /<?£  ),  which  can  be  useful  in  determining  the 

®BB  ®NB  A  A 

applicability  of  the  combiner  developed  above.  For  example,  if  0^  and  eBB 
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Figure  D-3.  Sensitivity  of  Combiner  to  Correlation  of  Estimates 
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are  expected  to  have  p  <  .  2,  then  as  long  as  10  log  <r£  /<r£  )  is  between  -5 

®BB  ®NB 

and  +5.  On  the  other  hand,  if  higher  correlations,  say  p  =  .  8  are  expected, 
then  the  combiner  will  be  useful  only  when  the  two  estimates  have  virtually  the 
same  variance  (and  even  then,  the  improvement  will  be  small). 
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APPENDIX  E 


COMPARISON  OF  A  SIN(X)/X  AND  QUADRATIC  INTERPOLATOR 

The  first  phase  of  the  tracker  study ^  ^developed  and  analyzed  an  interpolator 
for  extracting  a  delay  estimate  from  the  adaptive  filter  weights  when  the  inter¬ 
array  delay  falls  between  the  discrete  adaptive  filter  taps.  This  was  done  by 
forming  a  continuous  interpolated  weight  function. 


M-l 

h.(t)  =  £  W(m)  f(t-mT  )  (E-l) 


where  W(m)  is  the  m^  adaptive  filter  weight,  T  is  the  sample  interval,  and  f  is 

8 

an  interpolation  function.  In  [  1] ,  the  function  used  was 


„„  _  sin  2  ir  Bt 
f(t)  =  2*Bt" 


(E-2) 


The  delay  estimate  was  determined  by  locating  the  peak  of  the  interpolated 
function  (E-l).  The  analysis  of  the  performance  of  this  interpolator  followed 
reference  [4]  in  mapping  the  fluctuations  in  the  weight  vector,  W,  into  fluctua¬ 
tions  in  the  determination  of  the  peak  location. 

In  practice,  it  would  be  desirable  to  reduce  the  number  of  computations 
required  for  interpolation  relative  to  this  sin  x/x  interpolator.  This  is  particu¬ 
larly  true  in  light  of  the  three  primary  drawbacks  of  the  sin  x/x  approach: 

(1)  for  optimality,  the  input  processes  should  be  ideal  band-limited 
random  processes 

(2)  every  weight  value  is  used  in  the  interpolation,  even  though  some  may 


be  dominated  by  fluctuations 


(3)  an  iterative  search  routine  is  required  in  order  to  estimate  the  zero 


crossing  of  the  interpolated  derivative  of  the  weights 
A  quadratic  interpolator  overcomes  these  problems.  The  quadratic 
interpolator  simply  selects  the  largest  weight  and  the  two  adjacent  weights, 
and  makes  a  quadratic  fit  to  the  data  to  find  the  peak .  A  performance 
analysis  follows,  along  with  a  comparison  of  the  two  interpolation  approaches 
based  upon  numerical  evaluation  of  the  expressions  developed. 

Suppose  on  a  particular  iteration,  weight  Wj  is  the  largest  of  the  weights 
of  the  adaptive  filter,  and  Wj  and  WJ+^  are  the  two  weights  adjacent  to  WJ# 

The  quadratic  interpolator  will  fit  a  curve 

f(Z)  =  aQ  +  ax  Z  +  a2Z2  (E-3) 

such  that  it  passes  through  the  points  Z  =  WjV  lf  Z  =  Wj,  and  Z  =  WJ+1»  For 
such  a  curve,  the  coefficients  of  (E-3)  are  given  by 

(E-4) 

f  Kn  *  Vi  -  2Wj]  <*-*> 


(E-6) 


where  T  is  the  algorithm  sample  rate.  Then  the  parameter  Z  represents  the 

S 

continuous  value  of  delay  along  the  tapped  delay  line.  The  delay  estimate  is  the 


value  of  Z  corresponding  to  the  peak  of  (E-3).  To  locate  this  peak,  differentiate 
(E-3)  and  set  it  to  zero,  yielding  the  delay  estimate 


A 

T 


WJ-1  -WJ+1 
WJ-1  +  WJ+1  -2WJ 


JT 


(E-7) 


Note  that  r  is  a  ratio  of  correlated  random  variables  with  non-zero  mean, 
and  that  a  general  determination  of  its  density  requires  knowledge  of  the  densities 
of  the  numerator  and  denominator.  Since  the  weights  are  usually  assumed  to  be 
Gaussian,  the  temptation  is  to  assume  that  the  numerator  and  denominator  are 
Gaussian.  This  leads  to  an  infinite  variance,  however,  even  though  simulations 
of  the  quadratic  interpolator  indicate  good  performance.  In  fact,  the  densities 
are  not  Gaussian,  since  they  are  conditional  on  die  selection  of  the  weights  in  the 
vicinity  of  the  peak.  For  example,  the  density  of  the  numerator  is 

P  (Wj-l  +Wj+i  "2Wj  I  Wj  >  WK,  K  *  J)  (E-8) 


This  density  is  not  Gaussian,  and  is  not  easily  determined.  However,  when  it 
can  be  assumed  that  the  true  peak  of  the  weight  function  lies  between  Wj_j  and 
WJ+1,  the  method  of  [  4]  can  be  applied  to  determine  the  variance  of  r . , 
When  the  true  delay  between  the  half  array  outputs  is  in  the  vicinity  of  Wj, 
then  the  variance  of  r  is  given  by 


(E-9) 


with  Z.  the  true  value  of  delay.  Based  upon  (E-3)  through  (E-6),  the  numerator 
and  denominator,  evaluated  at  Z  ,  are  given  by 


dZ 


dm 


dZ 


E  [W.T-1]  +E  "  2E  tW.T] 


(E-10) 


and 


Var 


df(Z) 

2 

aw 

,  + ,,  (-L.)2 

dZ 

Z-Zt  2TS2 

1  “'V 

2  . 
r  11 

based  upon  [1].  Then 


(E-ll) 


Here,  ow  is  the  variance  of  the  weights,  assumed  constant  over  all  weights 


Var  (rj  = 


„  2t2 
°W  8 

1  *  12  ('t‘ )  2 
s 

2 

E^J-l 

)  +  E(WJ+1)  -  2E(Wj) 

(E-12) 


To  proceed  further,  statistics  are  needed,  in  particular,  the  input  correlation 
function.  These  are  available  from  Appendix  V  of  [1].  Assume  that 


s3(f) 


sa(f) 


Ps’  -Bs^  f  <  Bs 
0,  otherwise 

Po'  -  B8  S  S  BS 

0,  otherwise 


(E-13) 


(E-14) 


E-4 


Then  from  [  1  ] ,  using  the  frequency  domain  model,  the  statistics  of  the  FFT 
of  the  weight  vector,  where  H(k)  is  the  weight  of  the  k^1  FFT  bin  of  the- frequency 
domain  model,  are 

ZL 


E[H(k)l  * 


.2  it 


A 


P  +  P 


,  0  <  k  s  J-l 


,  J-l  <  k  s  M-J 


(E-15) 
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,  M-J  <  k  <  M-l 


and 


Var  [H(k)]  = 


(Ps+PN)(2^M(Ps.PN)r  *  0  S  K  *  J-l  or  M-J  <  k  *  M-l 


,  J-l  <  k  <  M-J 


The  mean  of  the  time  domain  weights  are 


(E-16) 


1  Hr1  J  If  m<krk2> 

EtWJ  =  £  E  [H(k)J  e  M  1  2 

m  k=0 

pa  •to  m-  (t;  -  m) 

■  i  Jfc  Y  a  ■ 


M  P_+P 


sin 


'(<- 


W'(H 


(E-17) 


E-5 


and 


E  |  W  |  2 
m 


Therefore, 

Var  [W  ] 
m 

Substituting  (E- 

Varfr]  * 


M-l 

E 

v° 


E  [  Hlkj)  ]  e 


1  Mm<kl 


2J  OMWV 

1?  <Vpn> 


17)  and  (E-19)  In  (E-12)  gives 


1  MpN(2Pa+PN)(Ps^PN)  r2 

2  P2  (2-J^M  (PS+PN))  q 


r 


r 


where 


< 


1 

Z 

m=-l 


2JMT 


sin 


sin 


3  ft  -I 


jr_ 

M 


ft-) 


T  COS 


(E-21) 


Recall  from  U1  that  the  variance  for  the  continuous  adaptive  tracker  and  the 
sin  x/x  interpolated  adaptive  tracker  can  be  written  in  the  same  form  as  (E-20) 
with  Kq  replaced  by  Kq  and  K^,  respectively,  where 


V  (Zt-mTs) 


cos  2  ir  B(Zt-mTg)  sin  2  tt  B(Zt-mTg) 
(Zt-mTs)  2  ir  B(Zt-mTg)2 


V  (0)  *  0 

2  sin  2  ir  B(Zt~mTs)  2  cos  2  ir  B(Zt-mTs) 
2irB(Zt-mTs)3  (Zt-mTs)2 

2  ir  B  sin  2  ir  B(Z.-mT  ) 

f"  (0)  =  - 

The  variance  of  the  peak  estimate  for  the  sin  x/x  and  quadratic  interpolators 
can  be  compared  by  an  examination  of  the  Kj  and  K^.  Therefore  it  is  only 
necessary  to  numerically  evaluate  the  expression  for  Kj,  (E-23),  and  compare 
it  to  the  expression  for  K^,  (E-21).  Whichever  term  is  larger  has  the  greater 
variance  and  consequently  worse  performance. 

The  expressions  for  Kj  and  were  programmed  on  a  computer  and 
numerically  evaluated.  A  sample  frequency  of  2400  Hertz  was  assumed  with 
a  16  tap  adaptive  filter.  A  value  of  2-10  was  used  for  p.  Since  the  values  of 
Kj  and  are  only  dependent  upon  signal-to-noise  ratio  to  the  extent  that  they 
are  evaluated  in  the  vicinity  of  the  true  signal  peak  of  the  weights,  an  arbitrary 
SNB  value  of  20  dB  was  used. 

The  value  of  Kj  was  evaluated  for  signal  bandwidths  of  300,  600,  900,  and 
1200  Hz  with  the  interpolator  bandwidth  parameter,  B,  matched  to  the  signal 


f’MZt-mTa)  = 


E- 


bandwidth  B  ,  in  all  cases.  It  was  also  found  that  K.  is  dependent  upon  the 
position  of  the  peak  within  the  filter,  so  signal  delays  of  1,  2,  4,  and  8  taps 
were  evaluated.  This  dependence  upon  peak  position  only  occurs  when  die  peak 
is  sufficiently  near  the  end  of  the  filter  that  part  of  the  peak  falls  outside  the 
filter.  Table  E-l  summarizes  the  values  of  Kj  for  these  parameter  selections. 

The  values  of  Kj  obtained  for  a  signal  bandwidth  of  300  Hertz  may  be  suspect. 
By  Equation  (10)  of  [l]  the  value  of  (J-l)  obtained  for  this  signal  bandwidth  is, 

(J-l)  =  [(300(1/2400)  (16/2)  J  »  1 

From  [  1]  the  expectation  of  the  weights  is: 

,  0  s  k  sJ-1 

,  J-l  <  k  <  M-J  (E-24) 

A(M-k) 

,  M-J  £k£  M-l 

Thus  for  (J-l)  equal  to  1  only  4  of  the  weights  have  nonzero  expectations.  This 
may  be  too  few  values  for  an  accurate  evaluation  of  Kj. 

Similarly,  was  also  evaluated  for  signal  bandwidths  of  300,  600,  900  and 
1200  Hertz.  In  the  derivation  of  the  expression  for  Kq  it  was  assumed  that  in 
the  worst  case  the  center  weight  used  for  the  quadratic  fit  would  be  located  1  tap 
away  from  the  true  peak  of  the  weights.  Therefore  the  value  of  A  in  (8)  will  vary 


TABLE  E-l.  THEORETICAL  VALUES  OF  Ki  FOR  VARYING 
SIGNAL  BANDWIDTHS  AND  DELAYS 


between  0.0  (for  the  center  weight  directly  on  the  peak)  and  1.0  (for  the  center 
weight  one  tap  away  from  the  peak).  Accordingly  bounds  on  the  value  of  were 
obtained  by  evaluating  E-21  with  A  equal  to  0.0  and  1.0.  Table  E-2  summar¬ 
izes  the  values  obtained  for  K  . 

<1 

Table  E-3  provides  a  comparison  of  the  performance  of  die  two  interpolators. 
By  forming  the  ratio  of  Kj/K^,  the  ratio  of  the  standard  deviations  of  the  two 
interpolated  estimates  can  be  directly  compared.  If  the  ratio  of  Kj/K^  is  greater 
than  1. 0  then  standard  deviation  of  the  sin  x/x  interpolated  estimate  is  greater 
than  that  of  the  quadratic  interpolated  estimate.  If  the  ratio  of  Kj/Kq  is  less  than 
1. 0,  then  the  sin  x/x  interpolator  has  the  smaller  standard  deviation. 


TABLE  E-2.  THEORETICAL  VALUES  OF  It,  FOR  VARYING 
SIGNAL  BANDWEDTHS  AND  £ 


— ——————— —1 

A«0.0 

A*  1.0 

Signal  bandwidth  *  1200  Hz 

3.5136  x  10"4 

4.9690  x  10"4 

Signal  bandwidth  ■  900  Hz 

3.9284  x  10"4 

5.5556  X  10"4 

Signal  bandwidth  3  600  Hz 

4.5361  x  10-4 

6.4150  x  10-4 

Signal  bandwidth  *  300  Hz 

5. 5556  x  10-4 

7.8567  X  10-4 

TABLE  E-3.  COMPARISON  OF  SIN  X/X  AND  QUADRATIC  INTERPOLATORS  -  THEORETICAL 
VALUES  OF  K./K  FOR  VARYING  SIGNAL  BANDWIDTHS  AND  DELAYS 


In  estimating  the  peak  of  the  weights  the  sin  x/x  interpolator  uses  an  iterative 

search  routine  to  estimate  the  zero  crossing  of  the  interpolated  derivative  of  the 

weights.  Since  the  quadratic  interpolator  requires  the  evaluation  of  only  a  single 

expression  in  order  to  estimate  the  peak  of  the  weights,  it  is  both  simpler  to 

program  and  requires  significantly  fewer  calculations  than  the  sin  x/x  interpolator. 

» 

Thus  it  appears  that  even  though  the  quadratic  interpolator  is  less  complicated 
than  the  sin  x/x  interpolator,  it  offers  advantages  in  both  performance  and  imple¬ 
mentation,  at  least  for  this  selection  of  parameters. 

As  a  check  on  this  conclusion  with  real  data  and  different  parameters, 

Section  4.0  presents  the  results  of  processing  acoustic  data  collected  at  sea 
using  a  line  array,  where  both  the  sin  (x)/x  and  quadratic  interpolators  were 
used.  The  adaptive  filter  used  to  process  this  taped  data  had  64  taps,  and  ft  was 
set  at  2  10 ,  2~12,  2~14,  and  2  16 .  The  input  data  bandwidth  was  480  Hz,  with  a 
sample  rate  of  1  KHz.  Using  the  results  of  this  section,  the  theoretical  per¬ 
formance  of  this  tracker  with  broadband  static  targets  has  been  computed  for 
this  range  of  p.  Figure  El  shows  the  predicted  performance  for  the  sin(x)/x 
interpolator,  while  Figures  E2  and  E3  show  the  predicted  performance  for  the 
quadratic  interpolator  when  A  =  0  and  A  *  1,  respectively.  As  above,  the 
A  =  0  and  A  -  1  case  bound  the  theoretical  bearing  variance  when  the  true  peak 
of  the  weights  is  within  one  tap  of  the  largest  weight.  For  this  set  of  parameter 
values. 


Kj/Kq  *  2.3  for  A  =  0 


and 


Kj/Kq  „  .  7  for  A  =  1 
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Figure  El.  Predicted  Standard  Deviation  of  Bearing  Estimate  for  Sea  Tape  Parameters  With 
Sinx/x  Interpolator  sample  frequency  *  lKHz.M  ■  64  Taps,  Bandwidth  *  480  Hz 
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Figure  E-2.  Predicted  Standard  Deviation  of  Bearing  Estimates  for  Sea  Tape  Parameters  With 
Quadratic  Interpolator  -  Sample  Frequency  *  lKHz,  M  *  64  Taps,  Bandwidth  *  480  Hz, 

A  ■  0 


i 


E-15 


22-01010 


Figure  E-3.  Predicted  Standard  Deviation  of  Bearing  Estimate  for  Sea  Tape  Parameters  with 
Quadratic  Interpolator  —  Sample  Frequency  -  1  KHz,  M  *  64  Taps,  Bandwidth  ■  480  Hz, 

A  -  1 


E-16 


I»'w. 


For  this  set  of  parameter  values,  the  sin  x/x  interpolator  can  provide  slightly 
better  performance  for  some  values  of  delay.  However,  the  computational 
efficiency  of  the  quadratic  interpolator  would  still  make  it  more  attractive. 
These  results  are  compared  to  those  achieved  with  the  sea  tapes  in  Section  4. 0. 
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APPENDIX  F 


CANCELLATION  OF  INTERFERENCE  PRIOR  TO  TRACKING 

There  are  basically  two  approaches  to  dealing  with  interference  in  the 
adaptive  tracking  system.  One  is  to  depend  upon  the  spatial  response  of  the 
split  arrays  and  upon  the  adaptive  properties  of  the  tracker  to  provide  relative 
immunity  to  the  interference.  The  other  approach  is  to  attempt  to  cancel  the 
interference  from  the  half  beam  inputs  to  the  adaptive  tracker.  It  was  shown 
during  phase  2  of  this  study  [  2]  that  the  primary  effect  of  an  interference  is  to 
introduce  a  bias  in  the  bearing  estimate,  and  that  this  bias  can  be  quite  severe 
when  the  signal -to-interference  ratio  is  low  at  the  split  array  outputs.  Under 
these  circumstances,  it  may  be  desirable  to  cancel  the  interference  prior  to 
tracking. 

The  adaptive  noise  canceller  [5]  is  often  used  to  eliminate  an  interference 
arriving  on  the  sidelobes  from  a  beam  output.  A  reference  sensor,  for  example 
an  omnidirectional  hydrophone  that  is  spatially  separated  from  the  array,  is  used 
as  a  reference  in  the  LMS  canceller  configuration  of  Figure  FI.  When  the  inter¬ 
ference  power  dominates,  the  adaptive  filter  will  spatially  reject  the  interference 
from  the  beam  output.  Depending  upon  the  spacing  between  the  reference  hydro¬ 
phone  and  array  and  upon  the  separation  of  the  target  and  interference,  some 
rejection  of  the  target  signal  will  occur,  but  this  will  usually  be  acceptable  in 
light  of  the  interference  rejection.  Figure  F2  shows  this  approach  applied  to  the 
adaptive  tracker.  A  reference  is  supplied  for  each  half  array  and  a  canceller 
implemented  in  each  half  array  output.  A  convenient  choice  for  a  reference  may 
be  a  single  hydrophone  located  in  the  opposite  half  array. 
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Figure  FI. 

LMS  Adaptive  Canceller 


Figure  F2. 

Adaptive  Cancellation  Applied  to  the  Adaptive  Tracker 


Inclusion  of  the  cancellers  of  Figure  F2  prior  to  the  tracker  will  add  the 
algorithm  noise  of  the  cancellers  to  beam  outputs.  In  a  well  designed  canceller, 
however,  this  noise  will  be  negligible  in  comparison  to  the  interference  power 
rejected,  and  will  therefore  be  tolerable.  A  more  important  possibility  is  that 
in  inserting  the  transfer  functions  necessary  to  cancel  the  interference,  the 
cancellers  may  modify  the  weight  vector  of  the  adaptive  filter  so  as  to  bias  the 
resulting  estimate  of  target  bearing.  For  broadband  targets,  this  can  be  analyzed 
by  determining  the  mean  weight  vector  of  the  adaptive  tracker  with  split  beam 
cancellation  and  observing  the  effect  of  the  cancellation  on  the  location  of  the 
peak.  With  narrowband  targets,  a  bias  will  be  reflected  in  the  phase  of  the  FFT 
bin  containing  the  signal. 

Analysis  of  this  effect  is  best  performed  in  the  frequency  domain  in  terms 
of  continuous  Wiener  filters.  In  steady  state,  the  adaptive  cancellers  converge 
in  the  mean  to  discrete  Wiener  filters,  which  will  approximate  the  continuous 
Wiener  filters  provided  sample  rate  and  filter  length  are  judiciously  chosen. 
Assuming  that  each  reference  is  located  in  the  opposite  half  array ,  referring 
to  Figure  F2,  the  inputs  to  the  cancellers  can  be  written 

*4(t)  -  /b<P,9a)s(t-Ts-p>  dp  +  /b(p,  9 j)  i  (t  -  Tj  -p)  dp 

+  <F-l> 

*r<t)  -  /b(p,  9a)  s  (t  +  t  8-p)  dp  +  / b(p,  9  j)  l  (t  +  Tj  -p)  dp 

+  V(t>  (F“2) 

yt(t)  "  *(«  t,  )  + 1  (t  -A  Tx  )  +  (F-J) 

9 

F-4 


yr<<>  -  +  i  )♦  V(t)  <F'*) 

where 

d  =  distance  from  split  array  phase  center  to  its  reference 
D  =  distance  between  split  array  phase  centers 
s(t)  =  signal 
i(t)  =  interference 

b(p,  9  )  =  impulse  response  of  array  to  plane  wave  source  at  bearing  0 
t  =  delay  between  phase  centers  for  the  signal 
t  j  =  delay  between  phase  centers  for  the  interference 

Here,  n  (t),  n  (t),  n  (t)  and  n  (t)  are  zero  mean  and  uncorrelated  with 
x  Jt  xr  yjt  yr 

each  other.  For  plane  wave  signal  and  interference. 
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and 


T  = 
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2C 


sin  6 
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(F^) 


where  6  g  and  0  j  are  bearings  of  the  target  and  interference,  respectively. 
Letting  (t)  and  Wf(t)  be  the  impulse  response  of  the  two  cancellers  after 
convergence  (i.e. ,  the  Wiener  filters),  then  the  inputs  to  the  adaptive  tracker 
are 

=■  xz(t)-  /Wt  (p )  y*  (t  -  p  )  d  p  (F-7) 
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and 


er(t)  =  *r(t)  -  /wr<P)  yr(t  -  P  )  dp 


(F-8) 


The  Wiener  filter  transfer  function  corresponding  to  the  adaptive  tracker 
weight  is  given  by 


Wj(io)  = 
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Vi 

S?  e 
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(F-9) 


where 

Sere4  («)  ■  cross  spectral  density  between  £r(t)  and  e^(t) 

and 

Se  e  (u)  »  spectral  density  of  e  (t) 

A>  * 


The  functions  Se  e  (<*>)  and  Sec  (w)  can  be  found  from  (F-7)  and  (F-8) 

*  X*  Ju 


using 


Sere4  («)  =  rj  E  [er(t)ejl(t+  t)JJ 


and 


Stih  (w)“  7  E  K(t>  V 


{E[et«>cl<,+T>]J 


where  7  [  J  denotes  the  Fourier  transform ,  then 

s  e  e  («)  -  sx  x  <w>  +  Wr  <w>  wi  <“>  Sv  v  <“> 
r  t  xrxi  r  Vi 


-  W  («)  S  (u)  -  W  (u)  S 
r  yri  ^  Vi 


(<•>) 


(F-10) 


F-6 


where 


=  y{E[Xr<t)  X*(t+T>]} 
sv  (»)  =  He  ha)y‘<*tT)]| 

*v  <•*-  yjs  ^r«)  >‘i<‘+t)]j 
8  (u>*  stJe  [*r(t>  y A  (t  +  r >J J 

with  W  (<jj)  and  W  (u)  the  transfer  functions  corresponding  to  W  (t)  and  W  (t). 

*  At  *  £ 

Also, 

sv.<">  *  %*,<“>♦  lw*(“,l2Vi"> 

-  2 Re  [wt«  Sv  M] 


(F-U) 


The  Wiener  filter  transfer  functions  associated  with  the  cancellers  are  given  by 


Sy  x  (<*>)  Sy  X  (“) 

W2(“>-S^T=)  - 

yiyi  Vr 


(F-12) 


Using  expressions  (F-l)  through  (F-4),  to  determine  the  required  spectral 
densities,  it  can  be  shown  that 


W*(“) 


B(u),  0g)S§(«)e  c  u  8  +  B(W,  0^)0  c  u  1 

PT<W> 


(F-13) 


This  can  be  simplified  to 


Sa(«)e  c  ( 

s  (u)  *  -S - - -  <  B(u),  6  J  Sjiw)  +  B(«,  9  )  <u>) 

er6e  PTZ(a.)  81  8  y 


-  B(o»,  0  jlSjiwJe 


-jfKxv T.) 


ST(w)e 


-j^I 


PT  («) 


B(w,  0  r)  Sg(u)  +  B(w,  0  j)  SN  (u>) 


-B(w,  0a)Sa(u)e 


(F-22) 

To  determine  the  adaptive  filter  weight,  Wf(w),  the  spectral  density  S  (w) 

EeEe 

must  also  be  evaluated  from  (F-ll).  Evaluation  of  the  required  spectral 
densities  yields 


P-10 


Se^(w)  =  0S»  2  Sg(a>)  +  |B(u,  6  j)|  2  Sj(w)  +  Sjj  (w)j  PT( 


-  |  B(u),  8a)  Sa(iif)e 


’3  c(X  "d)ts  l  ve  ,  "i'^(1‘D)Tal 


B(u,  epSjtuJe 


PT(“») 


|B(w,  0g)|  2  Sa(oj)  JS^oj)  +  SN  (u)J  +  |B(w,  0j)|2  Sj(u)  jSB(a))  +  ^  (oj)j 
+  Sjj  (u)  PT(u>)  -  2B(o),  0g)B(<u,  0  j)  Sa(cu)  S^u)  cos  [-£  (l  -  (ts  -  tJ 


PT(w) 


(F-23) 


Equations  (F-21)  and  (F-23)  can  be  used  to  write  a  general  expression  for  the 
Wiener  filter  transfer  function  corresponding  to  the  adaptive  filter  weight. 


'JI2ts 

s  (»)  e 

wtW-  ■  p-w- 


B(u,  ea)Sj(u»)  +B(0),  ©gJSjj  (u>) - B(w,  0j)  Sj(u})e 


Sj(u>)e 


-J  S2TI 


px<“) 


-ja(i-S)(x-T  )i 

B(w,  9  j) Sg(ta>)  +  B(u>,  0  j)  Sjj  (w)  -  B(<u,  0g)  Sg(<u)e  C  D  1  8  J 


(F-24) 


where  D(u)  is  the  numerator  of  (F-23) 

For  convenience,  let  Bs  ■  B<w,  Qa),  Bj  =  B(u,  0j),  andv  «  Sj,  (1-  1)  2  (^-Tj). 
Dropping  the  argument,  u  ,  define 


A*  IB  -B  e"jY|  =  Ib.-B  e  ^  I 
si  is 


=  tan 


-1 


Bj  sinY 


Bfl  -  Bj  cosTj 


and 


4>2  =  -tan 


-1 


B  sin^ 

8 

B,  -  B„  cosY 
IS  J 


In  this  notation,  (F-24)  can  be  written  as 


W.(u>) 


'  ¥.a*+s.sp_  Ib,I2+sisn  IbiI2  +  3n  p- 


S  -J  -2t 

_8  1  C  TS 

P  6 
T 


V"1*1 


2  o  1  Oi  a 

1 


p_.  PgA®  J  +®N  B[ 


+  ^T  e 
T 


(F-25) 


P-12 


If  the  interference  is  in  the  sidelobes,  then  Bg  »  B  ,  so  that  A  =» 


If  the  interference  dominates,  Sj  »  Sg  and  Sj  »  ,  then  P^,  as 


W{(a>) 


-j||2V2fl 


J  +s-  e 
SI 


J*2  ♦ 

Bs 


where,  from  above,  (Bj/Bg)  «  1. 


B  ,  and 

B 


(F-26) 
Sj.  and 

(F-27) 


P-18 


o 

Now,  assume  that  j  Bj  |  SN  «  SN  ,  since  involves  the  noise  from 

2  y  *  * 

the  entire  half  array  and  |  Bj  j  represents  a  sidelobe  response.  Then  for  high 
SNR  at  the  beamformer  output,  Bg  Sg  »  ,  (F-27)  reduces  to 


-j  [« 

W{(u)*se  lc 


2t8-2^ 


2Tj+  24.2 


(F-28) 


Note  that  if  the  cancellers  were  not  used  on  the  split  beams  this  would  have  been 


IB  I  2  S  -j  -  2t 
Wfa  (M)J  a  2  3  e  c  8  +e 

°  |B,I  S, 


-i  |  *Tl 


(F-29) 


On  the  other  hand,  for  small  SNR  at  the  beamformer  outputs,  B  S  «  SXT  , 

s  s  N 


Wf(u)%  S 


B.  I  2  S. 


N. 


“if— 
J  c 


fi2i 


SI 


SN' 


2, 


B. 


s 


S. 


s 


((F-30) 


while  the  transfer  function  in  the  absence  of  the  cancellers  would  be  as  in  (F-28). 
Therefore  the  overall  amplitude  of  the  weight  has  been  reduced  from  the  high 
SNR  case,  but  the  interference  has  still  been  suppressed  relative  to  the  signal. 

First,  consider  the  narrowband  tracker,  in  which  the  bearing  estimate  is 
extracted  from  the  phase  of  the  Fourier  transform  of  the  weight  sector  at  the 
signal  frequency.  From  [2],  it  can  be  seen  that  if  tracking  is  to  be  possible 
with  acceptable  accuracy  In  the  absence  of  the  interference,  then  BSS8/SN  »  1. 
Since  Sjj  includes  the  contribution  of  all  the  phones  in  the  array,  , 


F-14 


so  the  condition  Bg  Sg  »  Bj  is  met,  and  (F-28)  is  the  relevant  transfer 
function.  The  phase  at  a  frequency,  o> ,  is 


As  (Sj/Sg)  -*■  oo  ,  there  is  a  residual  bias  of  2^1  in  the  phase.  Note, 

however,  that  when  the  interference  is  in  the  sidelobes ,  B_  »  BT, 

3  1 


so  the  residual  bias  is  very  small.  It  can  be  made  arbitrarily  small  by  placing 
the  references  very  close  to  the  phase  centers,  so  (1  -  d/D)  «  0. 

Letting  ^  2Xg  -  2  ^  1  +  by  a  trigonometric  Identity, 


C^=  tan 


-1 


.ln[e2(Ti-T1)-2(  V  »2) 

r  f2)] 


(F-33) 


The  bias  in  the  phase,  B<J>,  is  then 

B,  ■  °i  -  2*i 


F-U 


Again,  when  Bg  »  BT,  it  can  be  seen  that 


«  c  (l-|)  2<Tj-ti)  »  *. 


(F-34) 


Substituting  this  in  (F-33)  for  Sj/Sg  »  1  gives 

B*  fe  ^  ,ln  c 

I  s 

s  j  Bt  j 

-5  2  -  _I  fi  -  “It  -  t  I 

2  ST  D  B  U  c  |Ts  TI! 

I  s 

Since  both  BT/B  and  S/ST  are  small,  the  bias  inserted  by  the  cancellers 
is  si 

is  small. 

In  the  broadband  adaptive  tracker,  a  bias  will  appear  as  a  time  shift  in  the 
weight  vector.  Consider  the  case  when  signal,  interference,  and  noise  all  have 
the  same  spectral  shape,  i.e., 


(F-35) 


Ss<U>  =PSSo(ai) 


Sj(w)  -  Pr  S0(u») 


(F-36) 


SN  <“>  "  PN  So(<J) 

X  X 


V” =  \  s°(“) 


tor  B,  Ps  «  P„ 


Wf(w)  = 


lB<*.  9a)rp3 


N 


2t, 


B(g>t8I)  PNy 
B(«,  9  s)  Ss 


(F-38) 


In  either  case,  the  interference  related  peak  will  have  been  sufficiently  suppressed 
to  have  negligible  effect  on  the  location  of  the  signal  peak.  Upon  inverse  trans¬ 
forming  to  obtain  the  mean  weight  vector  for  the  tracker,  using  (F-32), 
the  signal  peak  is  seen  to  occur  at 


2*s-2b;  I  (1-d>  2i*,-  ’il 


so  the  bias  in  delay  is,  substituting  from  (F-5)  and  (F-6) 


g 

Bt=  -2  r  "2  (1  ‘  |sin  ^  9I! 

s  C 

Again,  since  Bj  «  Bgf  the  bias  is  small,  and  can  be  made  arbitrarily  small  by 
placing  the  references  at  the  array  phase  centers. 
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size  on  bearing  lag  are  analyzed.  Adaptive  cancellation  of  interferences 
prior  to  the  tracker  is  shown  to  produce  negligible  bias  in  the  bearing 
estimate.  Bearing  estimator  structures  for  use  with  combined  broadband 
and  narrowband  targets  are  developed.  The  tracker  is  shown  to  track 
actual  sonar  targets  from  sea  tapes  with  variance  in  agreement  with 
theoretical  predictions.  A  single  tracker  structure  can  be  used  to  track 
multiple  targets  using  both  broadband  and  narrowband  energy. 
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